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When  studied  during  and  after  curing,  the  only  a  single  reaction  path 
was  found  during  curing,  but  the  reaction  path  was  highly  teirperature 
sensitive.  Physical  aging  is  found  to  be  an  increase  in  density  of 
amorphous  materials  that  must  occur  below  Tg  and  is  reversible  by 
heating  the  material  above  its  Tg.  The  aging  process  was  faster  if  it 
occurred  during  the  curing  reaction.  Fracture  energy  measurements  as  a 
function  of  surface-constrained  thin  film  thicknesses  show  that 
physical  constraints  on  the  plastically  deformed  crack  tip  may  reduce 
strength.  Creep  measurements  on  Viton  fluorohydrocarbon  elastomers 
crosslinked  to  different  degrees  confirm  that  looser  molecular  networks 
have  broader  viscoelastic  spectra  with  some  mechanisms  with  longer 
characteristic  times.  A  series  of  urethan  crosslinked  polybutadienes 
were  found  to  exhibit  an  exceptionally  low  lying  plateau  in  the  creep 
compliance  curve  which  had  not  previously  been  observed. 


TABLE  OF  CONTENTS 


Page 


OBJECTIVES .  1 

STATUS .  1 

CUMULATIVE  LIST  OF  PUBLICATIONS .  6 

THESES/DISSERTATIONS .  9 

INTERACTIONS .  11 

APPENDIX: 

I.  EPOXY  RESIN  STUDIES  OF  FRACTURE  IN  THIN  FILMS..  12 

II.  THE  EVOLUTION  OF  THE  VISCOELASTIC  RETARDATION 
SPECTRUM  DURING  THE  DEVELOPMENT  OF  AN  EPOXY 
RESIN  NETWORK .  2  6 

III.  EPOXY  RESINS  (DGEBA) ,  THE  CURING  AND  PHYSICAL 

AGING  PROCESS .  57 

IV.  VOLUME  DEPENDENT  RATE  PROCESSES  IN  AN 

EPOXY  RESIN .  88 

FINITE-ELEMENT  NUMERICAL  SIMULATION  OF  TEARING  AND 

FATIGUE  IN  POLYMERS  (Michael  R. Frank) .  109 


FINAL  TECHNICAL  REPORT 

"Fundamental  Studies  of  Time  Dependent  Response  and 
Fracture  of  Cross-Linked  Polymers" 

(a)  Objectives : 

i)  To  determine  the  relationships  of  polymer  structure, 
network  topology  and  other  microstructural  features  to  time- 
dependent,  small  deformation  properties  such  as  creep  and  volume 
relaxation  in  thermosetting  resins. 

ii)  To  determine  the  relationships  between  polymer  network 
structure  and  the  fracture  energy  of  thermosetting  resins. 

(b)  Status : 

i)  Epoxy  resins  based  on  diglycidyl  ether  of  bisphenol  A 
(DGEBA)  and  diamino  diphenyl  sulfone  (DDS)  have  been  prepared 
using  an  homologous  series  of  5  prepolymer  molecular  weights  and 
a  stoichiometric  ratio  of  reactants.  Resultant  networks  were 
characterized  as  a  graded  series  having  increasing  network  chain 
lengths  and  decreasing  Tg . 

ii)  The  series  of  Epoxy  resins  based  on  the  diglycidyl 
ether  of  bisphenol  A  has  been  extensively  investigated  during  and 
after  curing.  The  degree  of  cure  was  monitored  by  four  different 
methods:  1)  density  increase,  2)  heat  of  reaction  evolution,  3) 

glass  temperature,  Tg ,  increase,  and  4)  infrared  absorption  by 
the  glycidyl  ether  group.  The  epoxy  resins  studied  indicated 
that  only  a  single  reaction  path  was  followed 
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during  the  curing  but  the  kinetics  were,  of  course,  highly 
temperature  sensitive.  When  curing  temperatures  were  chosen 
that  were  lower  than  the  glass  temperature  of  the  fully  cured 
resins,  the  reaction  slows  down  when  the  Tg  increases  to 
temperatures  that  are  higher  than  the  curing  temperature,  Tc , 
but  the  curing  reaction  continues  more  or  less  indefinitely.  In 
about  a  month  the  Tg  can  become  50  degrees  higher  than  the 
curing  temperatures.  When  the  curing  is  slowly  progressing 
below  Tg  physical  aging  occurs  concurrently.  Physical  aging  is 
simply  an  increase  in  density  of  amorphous  materials  that  must 
occur  at  temperatures  below  Tg .  This  drift  toward  equilibrium 
density  is  accompanied  by  a  slowing  down  of  all  rate  processes 
occurring  in  the  material.  Physical  aging  and  all  of  its 
consequences  are  reversible?  i.e.,  its  effects  can  be  erased  by 
heating  the  material  above  its  Tg . 

It  was  determined  that  the  same  degree  of  aging  was 
achievable  after  the  curing  was  completed,  but  the  process  took 
less  time  if  the  aging  was  allowed  to  occur  during  the  curing 
reaction. 

The  curing  process  was  followed  not  only  by  determining  the 
extent  of  reaction  by  means  of  volume  shrinkage  and  by  changes 
in  infrared  absorption,  but  by  changes  in  the  viscosity,  gel 
fraction,  fictive  temperature  and  equilibrium  compliance. 

Reduced  creep  compliance  J(t)  curves  of  fully  cured  samples 
of  a  series  of  Epon  epoxy  resins,  where  the  molecular  weight  per 
cross-linked  unit  M*  varied  from  420  to  2870,  were  determined 
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and  their  retardation  spectra  were  calculated.  The  results 
showed  the  systematic  increase  in  the  number  of  long-time 
mechanisms  with  the  increase  in  M* . 

The  kinetics  of  physical  aging  as  monitored  by  the  decrease 
in  rates  at  several  temperatures  below  Tg  was  determined  for 
samples  with  four  different  crosslink  densities.  A  reduction 
procedure  was  devised  to  enable  comparison  of  the  results 
obtained  on  the  differing  samples. 

Different  stages  of  the  molecular  network  development  were 
prepared  and  from  the  measured  creep  compliance  curves  the 
evolution  of  the  viscoelastic  retardation  spectrum  that  occurs 
during  the  curing  of  an  epoxy  resin  was  determined.  It  was 
shown  how  the  short-time  behavior  at  Tg  is  independent  of  the 
material  and  the  stage  of  network  development;  but  the  long-time 
behavior  showed  a  dramatic  increases  in  retardation  mechanisms 
to  an  operationally  infinite  number  at  the  point  of  incipient 
gelation,  which  was  followed  by  an  equally  dramatic  decrease  in 
long-time  mechanisms. 

The  specific  volume  was  measured  for  one  of  the  Epon  epoxy 
resins  at  rates  of  cooling  which  varied  over  three  orders  of 
magnitude.  The  expected  decrease  in  Tg  with  decrease  in  rate  of 
cooling  was  observed.  The  ratios  of  these  rates  were  compared 
with  temperature  shift  factors  obtained  1)  from  the  kinetics  of 
volume  contraction  following  temperature  jumps  of  2.5°C,  2)  from 
torsional  creep  measurements,  and  3)  elongational  creep 
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measurements.  All  values  from  the  four  different  kinetic 
processes  fell  on  the  same  curve. 

iii)  Well-characterized  epoxy  networks  have  been  studied 
as  constrained  thin  films  in  mode  I  fracture  measurements. 
Fracture  energy  measurements  as  a  function  of  surface- 
constrained  thin  film  thicknesses  have  shown  that  physical 
constraints  on  the  plastically  deformed  crack  tip  may  well 
reduce  strength.  Also  thin  films  (<  0.5  mm)  may  be  expected  to 
show  the  effects  of  adherend  surfaces  as  reductions  in  the  state 
of  cure  (cross-link  density)  and  results  of  fracture  testing 
must  be  corrected  for  such  effects. 

iv)  Creep  measurements  on  three  Viton  f luorohydrocarbon 
elastomers  crossl inked  to  different  degrees  have  confirmed  that 
looser  molecular  networks  have  broader  viscoelastic  spectra  with 
some  mechanisms  with  longer  characteristic  times. 

A  corresponding  greater  spread  in  the  tear  energies  that 
was  expected  was  not  significant.  Tear  energy  data  as  a 
function  rate  were  reducible  with  the  temperature  shift  factors 
obtained  from  creep  compliance  measurements. 

However,  the  tear  results  do  appear  to  substantiate  the 
Mx1^2  proportionality  as  predicted  by  the  theory  of  Thomas  and 
Lake.  Mx  is  the  molecular  weight  per  crosslinked  unit. 

v)  A  series  of  urethane  crosslinked  polybutadienes  were 
found  to  exhibit  an  exceptionally  low  lying  plateau  in  the  creep 
compliance  curve  which  had  not  previously  been  observed.  Since 
the  presence  of  this  plateau  gives  rise  to  an  exceptionally 
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broad  viscoelastic  retardation  spectrum  higher  average  energy 
dissipation  is  seen  as  a  function  of  frequency.  Correspondingly 
broader  tear  energy  curves  are  also  observed.  Sustained  high 
losses  are  often  desirable  and  an  increasing  number  of  means  to 
create  such  losses  are  needed.  A  unique  strong  temperature 
dependence  was  also  found  to  be  exhibited  by  these  unusual 
elastomers. 

The  creep  and  tear  measurements  obtained  on  the  widely 
varied  elastomers  constitutes  a  unique  data  base  that  has 
enabled  us  to  correlate  the  linear  viscoelastic  behavior  with 
the  kinetics  of  failure.  An  equation  enabling  the  calculation 
of  the  rate  dependence  and  the  magnitude  of  the  tear  energy  from 
the  viscoelastic  retardation  function  has  been  found.  The 
latter  is  obtained  from  the  creep  compliance  curves.  Creep 
measurements  are  more  precise,  quicker,  and  more  convenient  to 
measure  than  the  tear  energy  curves.  In  addition,  the  under¬ 
standing  of  the  molecular  basis  for  creep  behavior  will 
therefore  be  transferable  to  rationalize  elastomeric  failure. 

vi)  The  integration  of  knowledge  gained  from  these 
studies  provides  strong  insights  into  the  origins  of  strength 
and  creep  properties  of  thermosetting  resins.  Also  the  cure 
path,  which  is  dependent  on  thermal  history  for  full  property 
development,  may  be  optimized. 
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(e)  Interactions 

Papers  presented: 

1.  "Fracture  of  Surface-Constrained  Epoxy  Resins," 

Monash  University,  Melbourne,  Australia,  June,  1989 

2.  "Relationship  between  Neat  Resin  and  Composite 
Properties,"  Keynote,  Asian-Pacif ica  Composites 
Conference,  Adelaide,  Australia,  June  1989 

3.  "Fracture  of  Surface  Constrained  Epoxy  Resins," 

Society  of  Plastics  Engineers  Conference,  Cleveland, 
Ohio,  October  1989 

4.  "Use  of  Parametric  Models  in  Designing  Polymeric 
Materials  to  Specifications,"  Society  of  Engineering 
Science,  Salt  Lake  City,  Utah,  September  1987 
"Correlation  of  Composite  and.  Bulk  Matrix  Resin 
Properties,"  Society  of  Plastics  Engineers,  Symposium  on 
High  Temperature  Polymers 
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I-  EPOXY  RESIN  STUDIES  OF  FRACTURE  IN  THIN  FILMS 


A.  Introduction 

The  fracture  properties  of  epoxy  networks  have  been  successfully  determined 
in  the  bulk[l-4],  as  adhesives[5-10],  and  in  fiber-reinforced  composites[ll-13] 
using  the  concepts  of  linear  elastic  fracture  mechanics.  However,  comparison  of 
the  data  from  the  various  types  of  test  specimens  indicates  that  the  bulk  network 
properties  do  not  translate  to  the  adhesive/composite  properties.  In  order  to 
better  understand  why  this  is  so,  a  logical  starting  point  is  the  realization  that 
there  is  a  major  difference  in  the  surface  to  volume  ratio  between  the  bulk  and 
the  adhesive/composite  type  specimens.  When  cast  as  adhesives  or  as  the  binding 
matrix  in  highly-filled,  fiber-reinforced  composites,  the  epoxy  networks  can  be 
envisioned  as  highly  surface-constrained  thin  films  which  have  very  high 
surface  to  volume  ratios.  Since  the  surface  layer  is  often  different  from  the  bulk, 
thin  film  properties  may  be  different  as  well. 

As  the  epoxy  networks  move  from  the  bulk  configuration  to  very  thin  films  the 
following  factors  may  become  increasingly  important: 

1)  Residual  stresses  arising  from  the  geometry  of  the  layer  and  bonding 
between  the  adherends/fibers  may  increase  so  that  the  resin  layer  is 
significantly  preloaded. 

2)  The  presence  of  the  adherends/fibers  may  chemically  or  physically 
affect  the  interfacial  (interphase)  layer.  Large  chemical  interactions 
may  cause  deviations  in  network  morphology  which  may  extend  well 
into  the  resin  layer. 

3)  The  high  modulus  adherends/fibers  may  constrain  the  plastically 
deformed  region  ahead  of  the  propagating  crack  causing  a  change  in  the 
amount  of  energy  dissipation  occurring  during  fracture. 

In  this  work  the  fracture  properties  of  the  epoxy  networks  cast  as  surface- 
constrained  thin  films  were  determined  and  compared  to  the  bulk  network 
fracture  properties.  Network  morphology  (average  chain  molecular  weight 
between  crosslinks,  Mc)  was  determined  for  epoxy  networks  cast  as  both 
constrained  (bonded)  and  unconstrained  (non-bonded)  thin  films. 

The  importance  of  knowing  Mc  becomes  obvious  based  on  previous  work  by 
LeMay(14]  which  showed  a  relationship  between  Mc  and  the  arrest  fracture 
energy,  Gica,  for  epoxy  networks  below  Tg  (see  Fig.  1).  Interestingly,  it  was  found 
that  the  arrest  fracture  energy  was  proportional  to  the  1/2  power  of  the 
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molecular  weight  between  crosslinks.  This  relationship  is  similar  to  that 
proposed  by  Lake  and  Thomas[15]  which  predicts  that  the  threshold  tearing 
energy  should  be  proportional  to  the  1/2  power  of  the  molecular  weight  between 
crosslinks  for  an  endlinked  network.  The  theory  requires  that  the  polymer 
chains  crossing  the  fracture  plane  have  enough  mobility  so  that  each  bond  along 
the  backbone  of  the  molecule  can  be  stressed  to  near  its  breaking  stress  before 
chain  breakage  occurs.  Althougth  the  epoxy  networks  are  well  below  Tg,  the 
occurrence  of  plastic  deformation  at  the  crack  tip  suggests  that  there  may  be 
sufficient  chain  mobility  occurring  for  the  theory  to  hold.  Thus,  variations  in  Mc 

caused  by  the  presence  of  adherends  or  reinforcing  fibers  in  thin  epoxy  layers 
could  be  responsible  for  discrepancies  in  fracture  energies  between  bulk  and 
adhesive/composite  configurations. 

Although  there  are  numerous  suggestions  of  a  'plastic  zone'  at  the  advancing 
crack  front  in  highly  crosslinked  epoxy  networks[l-3, 7-9,1 1-13,16,17], 
experimental  determination  of  the  size  and  shape  of  the  zone  is  extremely 
difficult,  if  not  impossible.  Yamini  and  Young[18]  measured  the  length  of  highly 
deformed  regions  of  material  on  the  fracture  surfaces  of  double  torsion  fracture 
specimens  at  positions  of  crack  initiation/arrest  in  an  epoxy  resin  and  found  that 
the  length  correlates  with  the  size  of  the  plastic  zone  predicted  by  the  classical 
Dugdale  model[19].  LeMay[14]  used  the  blunt  crack  analysis  model[10,20]  to 
calculate  the  critical  length  for  an  homologous  series  of  epoxy  resins  which 
differed  primarily  in  Mc  and  found  that  the  critical  length  was  equal  to  the 
plastic  zone  radius  predicted  by  the  classical  Irwin  model[21].  The  results  also 
showed  that  the  log  of  the  critical  length  varied  as  the  1/2  power  of  the  log  of  Mc 
(see  Fig.  2).  It  seems  reasonable  to  assume  that  the  plastic  zone  size  may  be 
dependent  upon  MCt  becoming  larger  as  Mc  increases. 

The  experimentally  determined  fracture  energy  for  epoxy  networks  is  found  to 
be  much  greater  than  the  energy  predicted  by  theoretical  calculations  based 
soley  on  bond  cleavage[22].  This  excess  energy  must  be  due  to  dissipative 
processes  occurring  in  the  material.  For  an  clastic,  glassy  polymer,  which 
fractures  at  small  strains,  it  seems  logical  to  assume  that  any  viscoelastic  and 
plastic  energy  dissipative  processes  occur  in  a  very  small  region  (plastic  zone) 
local  to  the  crack  tip.  If  the  development  of  this  zone  is  inhibited  in  any  way 
during  crack  growth  then  the  fracture  energy  may  be  different.  When  epoxy 
resins  are  used  as  adhesives  or  as  the  matrix  for  a  highly-filled  fiber-reinforced 
composite  the  layer  thickness  may  be  less  than  the  size  of  the  plastic  zone  (see 
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Fig.  3).  If  the  plastic  zone  size  (and  its  overall  volume)  is  diminished  by  the 
constraints,  the  fracture  energy  should  decrease. 

An  experimental  attempt  was  made  to  physically  constrain  the  plastic  zone  and 
measure  the  resulting  fracture  energy  by  using  a  modified  version  of  the 
’Outwater'  double  torsion  test[ 23-25J.  The  constrained  testpiece  is  similar  in  design 

to  one  suggested  by  Gent[26]  and  consist  of  a  thin  layer  of  epoxy  sandwiched 
between  two  aluminum  adherends  (see  Fig.  4).  Six  different  epoxy  networks 
which  varied  in  Mc  and  presumably  in  plastic  zone  size  were  used  in  this  work. 

The  thickness  of  the  epoxy  layer  was  systematically  decreased  from  approximately 
2.0  to  0.1  mm  and  the  resulting  fracture  energies  were  determined.  The  fracture 
energies  of  bulk  double  torsion  testpieces  were  also  measured  for  comparison. 

The  average  molecular  weight  between  crosslinks  was  measured  as  a  function 
of  film  thickness  for  both  constrained  and  unconstrained  epoxy  films. 

Constrained  (bonded)  films  were  cast  between  phosphoric  acid  anodized 
aluminum  sheets  which  were  subsequently  etched  away  completely  in 
hydrochloric  acid  after  the  network  had  cured.  Unconstrained  (nonbonded)  films 
were  cast  between  silicone  rubber  plattens  which  had  extremely  smooth 
'glasslike'  surfaces.  The  film  thicknesses  ranged  from  0.1  to  2.0  mm  for  the 
unconstrained  films  and  0.3  to  0.8  mm  for  the  constrained  films. 

B.  Experimental 
1.  Materials 

Stoichiometric  endlinking  of  difunctional  epoxy  resins  of  different  initial 
molecular  weights  with  a  tetrafunctional,  aromatic  amine  allowed  accurate 
control  of  Mc.  An  homologous  series  of  commercial  DGEBA  type  resins  (Shell 
Chemical  Company)  having  average  molecular  weights  ranging  from  370  to  4200 
g/mol  was  used  in  this  work.  The  curing  agent  used  was  4,4'-aminodiphenyl 
sulfone  (Aldrich  Chemical  Company).  The  resulting  networks  arc  identified,  in 

order  of  increasing  Mc,  by  the  trade  name  of  the  resin  from  which  they  were 
made:  828,  1001F,  1002F,  1004F,  1007F,  and  1009F.  The  higher  molecular  weight 
resins  were  synthesized  by  fusion  of  the  828  oligomer  and  thus  have  the  'F 
designation.  The  adherends  for  the  constrained  testpieces  were  the  aluminum 
alloy  6061-T6  (Alcoa). 

2.  Fracture  energy  measurements 

All  of  the  double  torsion  testpieces  were  fractured  at  25°C  on  an  Instron 
Universal  Testing  machine  at  a  crosshead  speed  of  0.5  mm/min.  The  compliance 
of  the  double  torsion  testpiece  varies  linearly  with  crack  length.  Thus,  the  stress 
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intensity  factor,  Kic,  is  proportional  to  the  applied  load  and  is  given  by  the 
following  relationship  (for  plane  strain) 

Kic  =  Pwm(3/wtnl3(l-v)]1/2  (for  w/2  >  6) 

where  P  is  the  load  at  fracture,  wm  is  the  torsional  moment  arm,  w  is  the  width  of 
the  testpiece,  tn  is  the  fracture  plane  thickness,  and  v  is  Poisson’s  ratio  (assumed 
to  be  0.35).  Details  of  the  mathematical  analysis  of  the  double  torsion  testpiece  can 
be  found  in  a  paper  by  Williams  and  Evans[25).  The  fracture  energies  of  the 
epoxy  networks  were  calculated  from  the  stress  intensity  factors  using  the  well- 
known  relationship 


Gic  =  Kic  2/e 

where  E  is  Young's  modulus  of  the  testpiece.  Young's  modulus  was  determined 
from  static  three-point  bend  measurements  based  on  simple  elastic  beam  theory. 

In  brittle  epoxy  networks,  cracks  often  propagate  intermittently  in  a 
'stick/slip*  manner  over  a  fairly  broad  range  of  test  temperatures  and  loading 
rates.  This  type  of  crack  growth  yields  a  load  vs.  time  (or  displacement)  curve 
which  is  sawtooth  shaped  (see  Fig.  5).  In  this  work,  all  the  epoxy  networks 
fractured  in  a  stick/slip  manner.  Therefore,  two  values  of  fracture  energy  were 
calculated  for  each  network.  Initiation  fracture  energies,  Gic>,  were  calculated 
from  the  average  value  of  the  crack  initiation  peak  loads  and,  likewise,  arrest 
fracture  energies,  Gica,  were  calculated  from  the  average  crack  arrest  values. 

3.  Mc  measurements 

The  average  molecular  weight  between  crosslinks  for  each  network  was 
determined  from  near  equilibrium  rubbery  tensile  modulus  measurements  made 
on  an  Instron  Universal  Testing  machine  equiped  with  an  environmental 
chamber  (see  Table  1).  The  moduli  were  determined  for  dogbone-shaped  epoxy 
films  ranging  in  thickness  from  0.1  to  2.0  mm.  The  stress-strain  measurements 
were  made  at  40°C  above  the  epoxy  glass  transition  temperatures.  Previous 

work[14]  has  shown  that  these  epoxies  behave  sufficiently  like  an  elastic  rubbery 
material  at  these  temperatures  and  reasonable  results  were  obtained  when  the 
kinetic  theory  of  rubber  elasticity  was  applied.  Me  was  calculated  for  each 
testpiece  using  the  following  relationship 
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Mg  =  3opRT/E 

where  0  is  the  front  factor,  p  is  the  density,  R  is  the  universal  gas  constant,  T  is 
the  absolute  temperature,  and  E  is  the  equilibrium  rubbery  tensile  modulus. 

C.  Summary  of  results 

1.  Mc  vs.  epoxy  film  thickness 

The  dependence  of  Mc  on  film  thickness  for  the  epoxy  networks  is  shown  in 
Fig.  6.  Data  was  not  obtained  for  the  828  or  1009F  networks  (the  828  testpieces 
were  too  brittle  to  handle  and  the  1009F  resin  was  too  viscous  to  be  cast  as  a  thin 
film).  No  variation  in  Mc  with  film  thickness  was  observed  for  the  1001F  network. 
However,  Mc  increases  as  the  epoxy  film  thickness  decreases  for  the  1002F,  1004F, 
and  the  1007F  networks.  The  deviation  in  Mc  becomes  greater  as  the  initial  resin 
molecular  weight  is  increased. 

The  variation  of  Mc  may  be  the  result  of  a  non-stoichiometric  reactant  ratio 
caused  by  a  loss  of  curing  agent  to  surface  contaminants  (e.g.  moisture).  Since 
the  surface  to  volume  ratio  is  extremely  large  for  very  thin  Films  as  compared  to 
bulk  castings,  it  is  conceivable  that  a  significant  percentage  of  curing  agent  may 
be  lost.  Small  quantities  of  resin  (approx.  100  g)  were  used  to  cast  the  thin  films. 

A  loss  of  as  little  as  330  mg  in  the  1007F  networks  would  result  in  a  10%  deviation 
from  a  stoichiometric  reactant  ratio  (epoxy  excess). 

Since  the  more  tightly  crosslinked  networks  require  larger  quantities  of 
curative  for  a  stoichiometric  reactant  ratio  for  the  same  amount  of  resin,  smaller 
deviations  (approx.  3%  for  the  100 IF  network)  would  be  expected  if  the  same 
amount  of  curing  agent  is  lost.  Previous  work[14]  has  shown  that  Mc  is  a  sensitive 
function  of  the  amine/epoxy  reactant  ratio  for  the  epoxy  networks  used  in  this 
investigation.  It  was  shown  that  deviations  of  20%  from  stoichiometry  (epoxy 
excess)  resulted  in  a  twofold  increase  in  Me  for  1002F  epoxy  networks. 

Thin  films  were  cast  between  anodized  aluminum  sheets  (6061-T6)  for  the  1002F 
networks.  No  difference  was  observed  in  Mc  as  compared  to  thin  films  cast 
between  silicone  rubber.  This  result  suggests  that  the  observed  variation  in  Mc 
with  film  thickness  is  not  due  to  residual  stresses  caused  by  bonding  and 
subsequent  thermal  shrinkage. 

2.  Gjca  vs.  epoxy  layer  thickness 

The  dependence  of  the  arrest  fracture  energies  on  epoxy  layer  thickness  (gap 
thickness)  obtained  from  the  constrained  double  torsion  testpieces  for  the  epoxy 
networks  is  shown  in  Fig.  7.  The  initiation  fracture  energies  are  not  shown.  The 
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initiation  values  mimic  the  behavior  of  the  arrest  values  but  are  greater  in 
magnitude.  Data  was  not  obtained  for  the  828  and  1009F  networks  (the  828 
testpieces  contained  numerous  cracks  caused  by  thermal  shrinkage  upon  cooling 
after  cure  and  the  1009F  resin  was  too  viscous  to  be  cast  between  the  aluminum 
adherends).  No  variation  in  Gica  was  observed  for  the  1001F,  1002F,  and  1004F 
networks.  However,  a  significant  decrease  in  Gica  was  observed  when  the  epoxy 
layer  thickness  was  reduced  to  less  than  0.5  mm  for  the  1007F  network.  It  is 
suggested  that  this  decrease  in  fracture  energy  may  be  due  to  physical 
constrainment  of  the  plastic  zone  at  the  crack  tip  by  the  aluminum  adherends. 

This  decrease  in  Gica  is  not  observed  in  the  more  tightly  crosslinked  networks 
because  their  plastic  zones  may  be  smaller  than  the  constraining  dimensions  of 
the  testpiece. 

Although  0.5  mm  is  much  greater  than  the  plastic  zone  dimensions  predicted 
for  the  bulk  networks  using  the  classical  plastic  zone  models,  it  should  be  noted 
that  the  stress  field  at  the  crack  tip  may  be  significantly  different  in  the 
constrained  configuration  than  in  the  bulk.  Strong  bonding  to  the  aluminum 

adherends  will  create  residual  stresses  in  the  epoxy  layer  which  may  increase  the 
size  of  the  plastic  zone.  In  the  bulk  double  torsion  testpieces  residual  stresses 
were  relieved  by  annealing.  The  constrained  double  torsion  testpieces  exhibited 
significant  residual  stresses  when  viewed  through  polarizers  set  at  extinction 
even  after  annealing. 

The  data  in  Fig.  7  have  been  corrected  for  the  change  in  arrest  fracture  energy 
caused  by  an  increase  in  Mc  with  decreasing  film  thickness.  The  data  in  Fig.  6 
were  fit  to  a  second  degree  polynomial  for  epoxy  film  thicknesses  less  than  1  mm. 
(There  is  no  known  significance  for  this  type  of  fit  but  it  did  provide  the  best  fit.) 
This  equation  was  used  to  calculate  Mc  for  the  appropriate  epoxy  layer 
thicknesses  obtained  in  the  constrained  testpieces.  The  arrest  fracture  energies 
were  corrected  using  the  half-power  relationship  found  by  LeMay[I4]  (see  Fig.  1). 
Although,  correcting  for  Mc  lowers  the  arrest  fracture  energy,  the  change  is 
about  40  -  50  J/m^  for  the  network  with  the  largest  Mc  variation  (1007F).  Thus,  it 
is  still  assumed  that  the  decrease  in  arrest  fracture  energy  for  epoxy  layer 
thicknesses  less  than  0.5  mm  in  the  1007F  network  is  due  primarily  to 
constrainment  of  the  plastic  zone. 

3.  Bulk  vs.  constrained  testpiece  fracture  energies 

The  initiation  and  arrest  fracture  energies  obtained  from  the  bulk  double 
torsion  testpieces  are  listed  in  Table  1.  The  values  are  the  averages  of  no  less  than 
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22  tcstpieccs  with  the  exception  of  1009F  (6  testpieces  averaged).  The  initiation 
and  arrest  fracture  energies  obtained  from  the  constrained  double  torsion 
testpiece  are  also  listed  in  the  table.  These  values  were  obtained  by  averaging  the 
fracture  energies  obtained  from  testpieces  in  which  the  epoxy  layer  thickness 
varied  from  1.0  to  2.0  mm.  The  initiation  fracture  energies  are  significantly 
higher  for  the  bulk  testpieces  as  compared  to  the  fracure  energies  obtained  from 
the  constrained  testpieces.  Good  agreement,  however,  is  found  between  the  arrest 
fracture  energies  from  the  two  types  of  testpieces.  This  suggests  that  the 
constraints  may  influence  the  crack  initiation  process  but  do  not  affect  crack 
arrest. 

4.  Unstable  fracture 

The  difference  between  the  the  initiation  and  arrest  fracture  energies  (Gic>  - 
Glca)  increases  as  Mc  increases  for  both  types  of  testpieces.  In  other  words,  the 
magnitude  of  the  'stick-slip'  behavior  increases  as  the  epoxy  networks  become 
tougher.  The  ratio  of  Gjci  to  Gica  is  shown  in  Table  1.  The  ratio  increases  as  Mc 
increases  for  the  bulk  double  torsion  testpiece  but  remains  fairly  constant  for  the 
constrained  testpiece.  Thus,  the  amount  of  'stick-slip'  behavior  is  found  to 
increase  dramatically  in  the  bulk  testpieces  as  Mc  increases  whereas  it  remains 
constant  in  proportion  to  molecular  weight  between  crosslinks  in  the  constrained 
testpiece. 
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Figure  3.  Crac):  tip  plastic  zones.  (a)  unconstrained  in  the 
bul,v  (b)  constrained  by  rigid  adherends  (c)  constrained  by 
hign  nodulus  fibers  in  a  composite. 


Figure  4. 
testpiece . 


Modified  version  of  the  ' Outvater '  double  torsion 


Figure  5.  Typical  load  vs.  tine  chart  traces  obtained  fron  the 
double  torsion  testpiece  for  epo;:y  fracture.  (a)  unstable  crack 
propagation  (b)  stable  crack  propagation. 
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Figure  7.  Plot  of  arrest  fracture  energies,  Gjca,  vs.  epoxy  layer  thickness 
obtained  from  the  constrained  double  torsion  testpieces  for  the  1001F,  1002F 
1004F,  and  1007F  netv/orks.  Dotted  lines  are  the  average  fracture  energies 
obtained  from  the  bulk  double  torsion  testpieces. 
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II .  THE  EVOLUTION  OF  THE  VISCOELASTIC  RETARDATION  SPECTRUM 
DURING  THE  DEVELOPMENT  OF  AN  EPOXY  RESIN  NETWORK 

INTRODUCTION 

Many  epoxy  resin  systems  are  formed  from  the  reaction  of  a  difunctional  epoxy  resin  with 
a  tetrafunctional  amine.  The  rate  of  reaction,  hence  the  curing,  depends  on  the  kind  of  diamine  used. 
Aromatic  amines  such  as  4,4/-methylene  dianiline  MDA  and  4,4'-diamino  diphenyl  sulfone  DDS 
react  more  slowly  with  diepoxides  than  aliphatic  diamines  and  yield  higher  glass  temperatures.  It  is 
well  known  that  three  reactions  are  usually  involved  in  the  curing  of  such  epoxy  resin  systems.1 

I  Reaction  of  the  primary  amine  hydrogens 

O  H  OH 

H2N-R-NH2+H2C  -  CH-R* - ^H2N-R-N-CH2-CH-R*- 


II  Reaction  of  the  secondary  amine  hydrogens 

H  OH  O  OH 

HjN-R-N-ChJc-R^+HjC^- \lH-R*-- j!,H2N-R-N-CH2-CH-R*- 

CH, 

I 

HC-OH 

I 

R* 


III 


Etherification 


OH  O 

I  /  \  k, 

H,N-R-N-CH,C-R*-+H,C  -  CH-R* — — 

I 

CH, 


H-C-OH 


R* 

I 


-R-N-CH,  CH-R*- 

I  I 

CH,  0-CH,-CH-R*- 

I  I 

HC-OH  OH 

I 

R* 

I 


If  the  rate  constants  k,  and  k2  differ,  the  intermediates  formed  on  the  way  to  a  full  cure  will 
vary.  This  variation  is  called  the  substitutional  effect.  With  k1  >  k2,  the  usual  case,  linear  chain 
growth  will  tend  to  precede  branching  with  subsequent  network  formation.  Nevertheless  with  exact 
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stoichiometric  concentrations  of  reactants  the  same  network  is  obtained  when  the  reaction  is 
complete.  The  etherification  reaction  is  usually  negligible  unless  an  excess  of  epoxide  group  is 
present  or  the  curing  is  carried  out  at  exceptionally  high  temperatures.  We  don’t  believe  we  have 
encountered  any  noticeable  etherification  at  curing  temperature  up  to  210*C  with  DDS  and  the 
diglycidyl  ether  of  bisphenol  A.2*3 

During  the  curing  of  an  epoxy  resin  the  viscoelastic  behavior  of  the  reactant  mixture  is 
expected  to  change  as  the  molecular  architecture  varies  during  the  reaction.  It  was  desired  to  follow 
the  evolution  of  the  retardation  spectrum  L  as  the  reaction  proceeded  from  the  mixture  of  reactants 
to  the  fully  completed  molecular  network. 

The  retardation  spectrum  can  be  determined  from  the  dynamic  shear  compliance  function 
J/(w)  and  ln(u)),  the  storage  and  loss  components  of  the  complex  dynamic  compliance  J*(w)  =  J7(w)  - 
iJ^w)  or  from  the  creep  and  the  recoverable  creep  compliances  J(t)  and  Jr(t).  All  of  the  compliances 
have  the  units  cm2/dyne  or  Pa'1.  The  pertinent  relations  are4: 

d  In  r  (I) 

-oo 

J//(W)  -  1  +  /  ^2r2  d  In  r  (2) 

'  -OO 

J(t)  *  Jr(t)  +  t/j;  ■  Jfl  +  /  L(l-e't/r)  d  In  r  +  t/rj  (3) 

-OO 

where  Jg  is  the  time  independent  glassy  compliance  which  in  principle  is  the  strain  per  unit  stress 
contributed  to  the  deformation  by  the  stretching  and  bending  of  bonds.  In  practice  it  usually  includes 
the  deformation  that  arises  from  very  local  molecular  motions  such  as  side  group  vibration  and 
rotations; 

t  is  the  time,  sec.; 

w  is  the  circular  frequency,  radians  sec.'1; 

t}  is  the  viscosity,  poise  or  Pa  sec.; 

and  r  is  the  retardation  time. 
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It  is  crucial  to  note  that  J"(w)  is  dominated  by  1  /ur)  at  low  frequencies  and  J(t)  by  l/t)  at  long  times 
in  the  terminal  zone  of  response.  Because  of  this  domination,  values  for  the  retardation  spectrum  at 
long  times  cannot  be  obtained  from  these  functions.  The  viscous  dissipation  of  energy  in  a  creep 
measurement  involves  the  mutual  displacement  of  molecules  past  one  another  resulting  in  a 
permanent  deformation;  whereas  the  time-dependent  recoverable  deformations  which  are 
characterized  by  L  reflect  the  orientation  imparted  to  the  molecules  during  the  approach  to  steady 
state  creep.  True  steady-state  creep  is  achieved  when  the  orientation  occurring  reaches  a  maximum 
limiting  value  under  constant  stress  conditions  and  not  by  an  experimentally  perceived  constant  rate 
of  creep,  as  is  sometimes  inferred.  The  time  dependent  motion  occurring  without  the  presence  of 
stress  during  recovery  is  thus  traceable  to  the  thermal  diffusion  of  the  viscoelastic  material  back  to 
a  state  of  maximum  entropy  with  isotropic  properties  and  randomly  disposed  molecules. 

It  was  not  possible  to  stop  the  reaction  at  various  stages  by  reducing  the  temperature  since  the 
viscoelastic  characterization  of  the  material  has  to  be  made  over  a  range  of  temperatures  where  the 
molecules  are  mobile  (above  Tg)  and  hence  usually  reactive.  Therefore  it  was  decided  to  complete 
the  reaction  of  all  epoxide  groups  at  various  stages  of  the  molecular  network  development. 
Subsequent  viscoelastic  characterization  by  means  of  creep  and  creep  recovery  measurements  would 
then  be  mrde.  The  different  stages  of  development  were  obtained  by  first  reacting  a  fraction  of  the 
epoxide  groups  with  the  monofunctional  amine,  methyl  aniline.  The  remaining  oxirane  rings  were 
then  reacted  with  a  stoichiometric  amount  of  DDS  in  the  frictionless  creep  apparatus  with  which  the 
creep  measurements  were  made. 

In  this  manner  our  structures  and  results  are  not  affected  by  any  substitutional  effect,  since 
the  secondary  amine  methyl  aniline  reacts  with  oxirane  rings  randomly  and  all  the  amine  hydrogens 
on  the  tetrafunctional  crosslinking  amine  DDS  are  reacted  before  the  creep  and  creep  recovery 
measurements  are  made.  Eight  different  concentrations  of  DDS  were  used  in  preparing  samples 
with  various  stages  of  network  development.  The  f ractions  of  the  equivalent  ratios  of  DDS  employed 
were  0.25,  0.35,  0.40,  0.45,  0.50,  0.60,  0.70  and  1.0. 
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EXPERIMENTAL 

Materials 

The  epoxy  resin  used  in  this  study  was  an  Epon  100 IF  produced  by  the  Shell  Corp.  It  was 
characterized  by  our  collaborators  LeMay,  Swetlin,  and  Kelley5  at  the  Polymer  Institute  of  the 
University  of  Akron.  Its  structure  is  given  by: 


ch3  an  ch  j 


Some  of  the  characterization  parameters  are: 

Mn  =  996  g/mole;  n  =  2.31;  and  it  fictive  temperature  Tf  =  31  C. 

The  tetrafunctional  crosslinking  diamine  used  was  4,4/  diamino  diphenyl  sulfone,  DDS,  (MW  *  248.3 
g/mole) 

0 

0 

with  a  melting  point  mp  range  of  17:°  -  178°C.  and  a  density  p(23°C)  =  1.38  g /cm3  and  the 
monofunctional  chain  stopping  amine  used  was  methyl  aniline  MA  (MW  =  107.2  g/mole)  with  a  mp\ 
=  -57*C  and  a  boiling  point  range  bp  ■  194  -  196°C.  The  reagent  grade  methyl  aniline  was  redistilled 
in  our  laboratory.  To  produce  a  series  of  partially  -  developed  epoxy  systems  representing  different 
stages  of  the  curing  process  which  reflects  different  stages  of  molecular  network  development,  we 
used,  as  curing  agents,  various  combinations  of  the  mono-  and  tetra-  functional  amines,  while 
maintaining  the  stoichiometry  between  the  oxirane  ring  concentration  and  the  amine  hydrogens. 
Sample  preparation 

As  an  example,  the  procedure  of  preparing  the  Epon  1001F/0.50  MA/0.50  DDS  sample  is 
given  below: 
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1)  A  sealed  flask  containing  the  mixture  of  8g  of  Epon  1001F  and  0.50  equivalent  ratio 
of  methylaniline,  was  heated  in  an  oil  bath  on  a  magnetic  stirrer  to  the  reaction 
temperature  of  143°C.  The  mixture  was  reacted  for  45  minutes  while  being  constantly 
stirred.  A  Fourier  Transformer  IR  (FTIR)  spectrum  showed  that  the  reaction  between 
Epon  100 IF  and  methylaniline  was  completed  in  less  than  45  minutes. 

2)  The  reaction  product  was  transferred  onto  a  Teflon  sheet,  cooled  to  the  room 
temperature,  and  was  ground  into  powder  with  a  porcelain  mortar  and  pestle. 

3)  0.50  equivalent  ratio  of  DDS  powder  was  added  to  a  measured  amount  of  the  reaction 
product  Epon  1001  F/0.50  MA.  The  well  -  mixed  powders  were  then  heated  to  160°C. 
When  the  DDS  completely  dissolved,  the  mixture  was  transferred  into  a  vacuum  oven 
at  110°C  and  was  degassed  until  it  became  a  bubble-free,  transparent  solution. 

4)  Portions  of  the  solution  were  transferred  into  DSC  sample  pans  which  were  placed  into 
the  sample  chamber  of  the  creep  apparatus.  These  samples  therefore  experienced  the 
same  thermal  h.jtoiy  as  the  material  whose  creep  behavior  was  determined.  The 
reaction  was  completed  in  the  creep  apparatus.  An  unchanging  equilibrium 
compliance  Je  measured  well  above  Tg  indicates  a  completed  reaction  as  well  or  better 
than  other  indicators. 

Instrumentation 

Fictive  temperatures  which  approximate  the  glass  temperature  Tf  2'3,6<7  were  measured  with 

'  *9 

a  Perkin  -  Elmer  Differential  Scanning  Calorimeter  DSC  II  (Main  Ave  -M.S.  1 2,  Norwalk,  CT  06856). 
Thermograms  were  obtained  with  a  heating  rate  of  10®/min  immediately  following  cooling  from 
above  Tg  at  a  rate  of  10°/min. 

Creep  and  recovery  measurements  which  yield  viscosities  and  recoverable  compliance  curves 
were  carried  out  with  a  frictionless  magnetic  bearing  torsional  creep  apparatus  (Time  -  Temperature 
Instruments,  P.O.  Box  40156,  Pittsburgh,  PA  15201)8.  Measurements  were  made  in  a  nitrogen 
atmosphere  after  drying  and  degassing  the  sample  in  vacuum  in  situ. 
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RESULTS 

Recoverable  Compliance  and  Viscosity  q£  Ih£  Neat  1001F  Resin. 

The  recoverable  shear  compliance  Jr(t),  (cm2/dyne  or  Pa'1)  and  the  shear  viscosity  coefficient 
J7,(poise,  dyne  sec/cm2  or  Pa  sec)  have  been  determined  for  the  neat  (unreacted)  Epon  100 IF  resin 
in  the  temperature  range  30°  to  77°C.  Jp(t)  curves  were  obtained  at  9  temperatures  from  30.2°  to 
53.5°C;  see  Fig.  1.  A  dramatic  10  fold  decrease  of  the  steady  state  recoverable  compliance  Je°  is  seen 
to  occur  with  a  20°C  decrease  in  temperature.  This  strong  temperature  dependence  is  shown  in  Fig.2. 

Consequently  the  log  Jr(t)  curves  cannot  be  ;  educed  to  a  common  curve.  In  Fig.  3  the  log 
Jr(t)  curves  are  shifted  along  the  logarithmic  time  scale  to  a  reference  temperature  T0  chosen  to  be 
30.2°C  by  using  the  logarithm  of  the  viscosity  ratio  r?(T0)/t?(T).  v(J0)  and  r?(T)  are  the  viscosities  at 
the  reference  temperature  and  at  the  other  temperature  of  measurement,  respectively.  This  plot 
attempts  to  compensate  for  the  difference  molecular  mobilities  at  the  different  temperature  and 
clearly  depicts  how  the  curves  fan  out.  This  thermorheological  complexity  persists  even  if  the  data 
are  reduced  to  a  common  Je°.  Reduced  recoverable  compliances  J*(t)  were  calculated  with  the 
expression 

JRr(t)  =  M(Jr(t)  -  Jg)  +  Jg  (4) 

where  Jg  =  the  glassy  compliance  =  8.0  x  10'11  cm2/dyne 

M  =  Jc°<To>  -  J9  <T0>  (5) 

Je°(T)  -  Jg  (T) 

T0  =  the  chosen  reference  temperature 

T  =  the  temperature  of  measurement. 

The  glassy  compliance  must  be  left  out  of  the  reduction  because  it  results  from  mechanisms 
found  at  shorter  times  that  have  a  much  smaller  temperature  sensitivity.  Jg  was  assumed  to  be 
temperature  independent  over  the  23°C  involved.  The  process  was  carried  out  with  the  result  shown 
in  Fig.  4  that  the  resulting  curves  were  not  superposable  because  the  higher  temperature  logarithmic 
curves  rise  more  abruptly  toward  log  Je°. 

Sixteen  viscosities  were  measured  between  30.2°  and  76.5  where  rj  fell  from  1 .32  X  1012  poise 
down  to  5.28  X  103  poise.  The  data  which  varied  over  the  eight  orders  of  magnitude  are  well 
described  by  the  Vogel,9  Fulcher,10  Tamman  and  Hesse,11  VFTH,  equation 
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(6) 

T-T„ 

rj  =  Ae 

or  log  t)  =  log  A  +  C/2.303  =  -2.36  +  470 
T  -  T.  T-0oC 

where  A,  C,  and  T.  are  characterizing  parameters.  The  latter  two  are  usually  interpreted  in  terms 
of  a  free  volume  model.4 
Curing  Simulation 

In  the  early  stages  of  curing  the  average  molecular  weight  increases  and  branching  appears. 
In  the  case  of  an  epoxy  resin  displaying  an  extreme  substitution  effect  linear  molecules  would  grow 
to  appreciable  lengths  and  hence  high  molecular  weights  before  branch  points  would  appear.  In  this 
study  branching  should  occur  simultaneously  with  molecular  weight  increase,  since  oxirane  rings  are 
randomly  reacted  first  with  the  secondary  amine  MA  and  all  of  the  DDS  hydrogens  are  reacted  before 
our  measurements  begin.  When  only  25%  of  the  oxirane  rings  are  reacted  with  DDS  crosslinking 
molecules  it  can  be  seen  in  Fig.  5  that  Je°  has  increased  by  about  30  fold  reflecting  the  increase  in 
molecular  weight.  Short  of  the  point  of  incipient  gelation  the  curing  system  is  a  viscoelastic  liquid 
whose  retardation  spectrum  in  the  terminal  zone  is  determined  by  the  molecular  weight;  more  so  by 
its  distribution,  and  by  long  chain  branching. 

It  should  be  noted  that  in  Fig.  5  for  the  system  with  25%  crosslinking  DDS,  and  for  all  higher 
concentrations,  the  thermorheological  complexity  exhibited  by  the  neat  resin  has  disappeared;  i.e.  all 
of  the  cured  systems  yielded  Jr(t)  curves  that  were  reducible  to  a  common  curve  within  the 
experimental  scatter  of  the  data.  Although  we  have  encountered  a  number  of  viscoelastic  liquids 
which  show  deviations  from  the  time-temperature  equivalence  principle4,  no  polymer  sample  with 
a  molecular  network  has  been  observed  by  us  to  show  such  deviations.  The  reducibility  of  all  of  the 
data  beyond  that  of  the  neat  resin  makes  our  presentation  simpler  and  clearer. 

At  35%  DDS  a  modest  intermediate  plateau  is  present  in  the  reduced  recoverable  compliance 
curve  as  seen  in  Fig.  6.  This  is  strongly  suggestive  of  an  entanglement  plateau.  If  a  bifunctional 
analogue  of  DDS  were  synthesized  this  conjecture  could  be  tested. 


34 


It  was  expected  that  the  terminal  value  for  Jp(t)  would  increase  to  infinity  at  the  point  of 
incipient  gelation  where  an  effectively  infinite  molecular  weight  was  reached.  Indeed,  when  an 
equivalent  ratio  of  0.45  DDS  was  used  in  the  curing,  the  reduced  creep  compliance  Jp(t)  =  J(t) 
[Tp/T increases  nine  orders  of  magnitude  from  the  glassy  level  before  indicating  a  tendency  to 
level  off  as  seen  in  Fig.  7.  Which  clearly  shows  it  is  just  beyond  the  gel  point.  The  reduced 
compliance  Jp(t)  is  used  for  temperature  shifting  of  crosslinked  systems  only.  It  has  been  shown  to 
be  inappropriate  for  uncrosslinked  polymers.12 

The  complete  effect  of  the  curing  of  this  epoxy  resin  is  depicted  in  Fig.  8  where  the  time 
dependent  recoverable  compliance,  J(t)  -  Jg  0,  is  shown  for  all  of  the  materials  cured  with  different 
reaction  ratios  of  amines  as  a  function  of  the  reduced  time.*  Jg  0  is  the  glassy  compliance  obtained 
by  extrapolating  the  short  time  data  to  zero  time  with  an  Andrade  plot. 

J(t)  =  J9,0  +  ^  +  */*  W 

where  J(t)  plotted  as  a  function  of  the  cube  root  of  time,  yields  the  glassy  compliance  from  the  zero 
time  intercept.  This  functional  form  is  normally  found  at  compliances  below  10'9cm2/dyne.  The 
time-dependent  recoverable  creep  compliance  that  is  common  to  all  of  the  systems  studied  is  most 
clearly  seen  when  the  glassy  compliance  is  not  included  in  the  comparison.  In  general  the  Jg  values 
should  be  analyzed  separately.  At  the  lower  degrees  of  crosslinking  (0.25  through  0.40  DDS)  the 
molecular  network  is  not  developed  and  hence  the  fully  reacted  materials  are  still  viscoelastic  liquids. 
The  dashed  lines  shown  in  Fig.  8  represent  the  recoverable  compliance  Jp(t);  whereas  at  the  greater 
degrees  of  branching  the  molecular  network  is  intact  and  the  solid  lines  represent  the  creep 
compliance  J(t).  Of  course,  for  a  crosslinked  material  J(t)  =  Jr(t).  The  reference  temperatures  T0 
used  in  the  preparation  of  Fig.  8  were  chosen  to  superpose  the  short  time  portions  of  the  J(t)  -  Jg 
curves.  It  can  be  seen  in  Table  I  that  the  T0s  required  are  on  the  average  two  degrees  higher  than 
the  measured  T#  temperatures.  We  consider  this  difference  to  be  within  the  uncertainty  of  our 
knowledge  of  corresponding  glass  temperatures.  The  evolution  of  the  recoverable  compliance  appears 


J  (t)-J  0  is  shown  for  the  viscoelastic  solids  with  molecular  networks  (solid  lines  in  Fig.  8) 
and  J  (t)  -  J.  _  is  shown  for  the  viscoelastic  liquids  (dashed  lines);  i.e.  specimens  with  no  gel 
fraction.  9,0 
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to  be  reasonable  and  straightforward.  Je°  increases  with  molecular  weight  (0.25  DDS).  A  rubberlike 
plateau  appears  with  further  increase  in  molecular  weight  and  branching  (0.35  DDS).  With  further 
branching  and  a  presumed  broadening  of  the  molecular  weight  distribution,  the  response  broadens 
on  the  time  scale  and  the  entanglement  plateau  disappears.  At  the  point  of  incipient  gelation  the 
macroscopic  tenuous  molecular  network  is  infinitely  compliant  (0.45  DDS).  As  the  network  further 
develops  and  becomes  tighter  the  enormously  long-time  mechanisms  disappear  due  to  the  intervention 
of  the  branchpoints  (0.50  DDS).  This  process  of  eliminating  long-time  mechanisms  which  involves 
cooperative  motions  over  long  distances  with  the  concomitant  decrease  of  the  equilibrium  compliance 
(0.60,  0.70,  and  1.0  DDS)  will  be  clearly  depicted  below  by  the  narrowing  of  the  retardation 
spectrum.  In  spite  of  enormous  changes  observed  in  the  long  time  deformations  with  the 
development  of  the  tight  molecular  network  the  short  time  deformation  appears  to  be  unchanged. 
Of  course,  the  position  of  the  compliance  curves  on  the  time  scale  is  dictated  by  the  local  free  volume 
which  decreases  with  the  elimination  of  dangling  chain  ends  and  the  increase  of  the  concentration 
of  tetrafunctional  branch  points.  This  decrease  in  free  volume  with  further  network  development 
is  the  reason  for  the  ascending  Tf  .  The  coincidence  in  form  of  all  of  the  recoverable  response  at 
short  times  reflects  the  fact  that  locally  the  polymer  chain  segments  are  the  same  and  the  short  time 
contributions  to  the  deformation  arise  only  from  the  local  segmental  chain  motions.  This  coincidence 
at  very  short  times,  may  in  addition,  just  reflect  contributions  to  the  deformation  which  are  universal 
for  all  glassy  materials. 

The  retardation  spectrum  obtained 13,14  for  the  earliest  stage  of  pre-gel  network-development 
from  the  data  shown  in  Fig.  5  (0.25  DDS)  is  presented  in  Fig.  9.  A  short  time  Andrade  region  with 
a  log-log  slope  of  1/3  is  followed  by  a  slightly  skewed  peak.  This  L(r)  is  typical  for  a  linear  polymer 
with  a  molecular  weight  which  is  less  than  Me,  the  molecular  weight  per  entangled  unit15,16. 

The  plateau  seen  in  Fig.  6,  which  is  presumed  to  reflect  the  presence  of  an  entanglement 
network,  gives  rise  to  the  development  of  the  usual  two  peaks  in  L(r)4,  as  seen  in  Fig.  10.  In  this 
case  the  first  peak  appears  as  a  strong  shoulder  preceding  the  rather  high  long  time  terminal  peak. 
This  shoulder  is  largely  filled-in  as  the  molecular  weight  distribution  is  presumed  to  broaden  with 
increased  branching;  (0.40  DDS).  This  can  be  seen  in  the  comparison  plot  of  all  of  the  retardation 
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spectra  obtained;  shown  in  Fig.  11.  The  divergence  of  L(r)  toward  infinity  at  long  times,  near  the 
gel  point,  is  clearly  seen  for  the  material  cured  with  45%  DDS.  L(r)  increases  ten  orders  of 
magnitude  from  the  glassy  level  of  response  up  to  the  value  at  long  times  where  the  material  is 
enormously  compliant  (ca  0.1  cm2/dyne). 

Both  comparison  plots  of  J(t)  and  L(r),  Figs.  8  and  11,  clearly  indicate  from  the  regions  of 
common  response  that  the  curves  for  reduced  times  less  than  109  seconds  show  characteristic 
curvature.  Only  at  longer  times  does  there  appear  to  be  exponential  behavior  which  is  indicated  by 
the  straight  line  segments  extending  for  approximately  5  measured  decades.  The  lcng-time  straight- 
line  slopes  are  0.75  for  d  log  Jr(t)/d  log  (t/at)  and  0.88  for  d  log  L(r)/d  log  rr. 

Temperature  Dependence 

The  temperature  dependence  of  the  shift  factors  aT  used  in  obtaining  all  of  the  reduced  creep 
and  compliance  curves  presented  in  Fig.  8  are  shown  in  Fig.  12  where  the  data  are  linearized 
according  to  variation  of  the  VFTH  equation 

aT  =  A'e  T'T«  (8) 

where  the  constants  A  and  C  are  1.231  x  10"16,  1850,  respectively  and  TB  =  Tf  g  -  50.5°C. 

This  equation  is  equivalent  to  the  Williams,  Landel  and  Ferry  equation17  and  therefore  is  an 
expression  of  the  free  volume  determination  of  rate  phenomena.  The  fit  of  the  data  was 
accomplished  in  the  manner  used  to  represent  the  aT  dependence  for  polystyrenes  of  widely  varying 
molecular  weights.  The  resin  cured  with  0.35  equivalent  ratio  of  DDS  was  chosen  as  the  reference 
material.  Its  optimum  T  was  50.5°C  below  its  measured  Tf  the  best  estimate  of  the  glass 
temperature  obtained  from  heating  experiments  with  a  differential  scanning  calorimeter  DDS2,3, 18, 19. 
The  heating  at  20°C/min  immediately  followed  cooling  from  above  Tg  at  the  same  rate*.  Its  Tf  was 
chosen  as  the  reference  Temperature  T0  where  aT  =  1.00.  Corresponding  Tms  were  chosen  for  all  of 
the  other  materials  so  that  the  same  constant  C  was  obtained.  It  was  then  assumed  that  the 
corresponding  "Tf  g"s  were  all  50.5'C  above  the  obtained  T^s.  With  these  "Tf  g"s  chosen  to  be  T0  as 


Therefore  little  physical  aging  took  place  while  the  material  was  below  its  T  . 
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the  value  for  A  *  1.231  X  10'16  obtains.  The  values  of  the  various  characterizing  temperatures  used 
in  this  and  the  compliance  analysis  presented  above  are  shown  for  comparison  with  the  DSC  Tf 
values  in  Table  I. 

The  viscosities  of  the  viscoelastic  liquids  0.25,  0.35,  and  0.40  DDS)  at  the  reference 
temperature  T0  for  the  reduction  of  the  compliance  curves  are  given  in  Table  II.  The  temperature 
dependence  of  the  viscosities  of  these  liquids  were  found  to  be  the  close  to  that  of  the  recoverable 
compliance  curves  over  the  following  temperature  ranges:  for  0.25  DDS,  65°  to  96°C;  for  0.35  DDS, 
70  to  122*C;  and  for  0.40  DDS,  79*  to  158“C.  Viscosities  can  be  calculated  with  Eq.  (6)  wh^e  the 
VFTH  parameters  are  given  in  Table  II.  All  of  the  values  for  C  are  the  same  as  that  from  the 
recoverable  compliance;  see  Fig.  12.  For  0.25  DDS  the  T„  of  15.8°C  is  the  same  as  that  for  Tp(t). 
For  0.35  and  0.40  DDS  the  Ta  =  10°C  are  a  few  degrees  lower  than  that  used  for  Jp(t)  indicating  a 
slightly  smaller  sensitivity  to  temperature  for  rj. 

It  is  of  interest  to  note  that  the  viscosities  determined  on  the  neat  1001°F  resin  were  best  fitted 
to  a  VFTH  plot  with  T«,  =  0°C  so  that  A  =  (Tf  g  -  TJ  *  31*C  and  C  =  1080.  These  values  for  C  and 
A  differ  substantially  from  the  values  of  1850  and  50.5°  determined  from  the  fitting  of  all  of  the  data 
for  Fig.  12.  However,  since  all  of  the  cured  systems  showed  no  measurable  change  in  Je  or  Je°  with 
temperature  as  did  the  neat  resin,  it  was  conjectured  that  the  temperature  dependence  of  the 
characteristic  time  constant  rc  ■  ij Je°,  which  is  a  pure  rate  parameter  with  seconds  for  units  (not  dyne 
sec/cm2,  which  are  the  units  of  the  viscosity  coefficient).  A  plot  of  log  f?Je°  versus  (T-TJ'1  with 
T„  =  -20°C  yielded  a  good  fit  to  all  but  the  lowest  temperature  points.  From  this  rc  fit  C  =  1890  and 
A  =  51*C.  These  values  are  certainly  within  experimental  error  of  the  values  obtained  from  the  shift 
factors  in  Fig.  12.  The  more  fundamental  nature  of  rc  as  a  rate  parameter,  rather  than  rj,  has  been 
noted  previously  under  similar  circumstances  where  Je°  of  a  low  molecular  weight  polystyrene  was 
found  to  be  strongly  dependent  on  temperature20. 

Real  Time  Curing  Bevond  the  Gel  Point 

With  the  stoichiometric  amount  of  crosslinker  DDS  **  1.0  used  in  the  curing  of  Epon  100 IF 
it  was  possible  to  periodically  obtain  J(t)  curves  from  about  one  second  out  to  several  hundred 
seconds.  In  all  of  the  eight  measurements  shown  in  Fig.  13  the  compliance  obtained  at  creep  times 
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greater  than  one  hundred  seconds  were  close  to  the  equilibrium  values.  The  approach  toward  zero 
slope  is  the  evidence  for  this  claim.  The  measured  J(t)  curves  are  indicated  by  the  solid  line  curves. 
The  long-dashed  lines  were  deduced  in  the  following  manner.  The  Tf  values  shown  in  Fig.  13  were 
obtained  from  the  relationship  established  previously;  see  Fig.  9  in  reference  2.  The  common  short- 
time  compliance  curve  at  Tf  was  obtained  from  Fig.  8  and  a  glassy  compliance  of  1.0  X  10‘10 
cm2/dyne.  Temperature  shift  factors  for  the  J(t)  curves  at  each  stage  of  cure  were  then  obtained 
from  the  common  temperature  dependence  shown  in  Fig.  10.  With  these  log  aT  values  each  curve 
wa:  shifted  from  its  Tf  to  the  temperature  of  cure,  142.8°C.  Where  the  deduced  curves  were  close 
to  the  short  time  measured  points  (at  22.4  and  120.6  hrs.  of  curing)  slight  addition  empirical  shifts 
were  allowed  to  obtain  reasonable  extrapolations  to  short  times.  These  are  represented  by  the  short- 
dashed  lines. 

To  make  possible  a  comparison  between  the  results  obtained  on  the  fully  reacted  but  partially 
developed  networks  and  the  reacting  and  developing  network,  the  four  extrapolated  curves  from  Fig. 
13  were  reduced  as  shown  in  Fig.  14.  The  curve  corresponding  to  22.4  hours  of  curing  was  chosen 
as  the  reference  curve  at  the  temperature  of  the  curing.  The  curve  with  the  Tf  g  of  96.8°C  represents 
a  degree  of  cure  of  62%;  T,  a  of  107.2°C  -  72%;  125.0°  -  88%;  and  135.1°  -  98%.  Examination  of 
Figs.  8  and  14  together  shows  that  the  extent  of  the  time  dependence  and  the  approach  toward 
equilibrium  compliances  are  comparable.  This  indicates  that  beyond  the  gel  point  the  networks 
formed  are  similar.  The  reduced  curves  in  Fig.  14  are  shifted  to  shorten  times  relative  to  those  in 
Fig.  8  by  about  four  decades  of  time.  This  is  consistent  with  the  fact  that  the  reference  temperature 
in  Fig.  14  is  approximately  68°  above  the  T„  of  the  reference  material;  see  Fig.  12,  (T-TJ'1-  14.7 
X  IO  3. 


There  has  been  considerable  interest  recently  in  the  viscoelastic  behavior  of  the  nascent 
molecular  network  that  forms  at  the  point  of  incipient  gelation.  The  behavior  at  high  compliance 
(beyond  rubbery  levels)  and  long  times  has  been  exclusively  examined.  Hopes  for  a  universal 
behavior  have  been  pursued(21)  and  percolation  theory(22'24>  suggests  material  independent 
exponential  behaviors  of  the  viscoelastic  functions.  Since  most  investigators  show  a  proclivity  to  only 
examine  and  consider  modulus  functions,  the  readily  interpretable  compliance  curves  and  retardation 
spectra,  that  have  been  obtained  in  this  investigation,  have  been  transformed  to  modulus  curves  and 
relaxation  spectra<25).  For  the  transformations  made  here  no  approximation  relations  were  used. 
Only  iterative  computer  calculations  were  used  for  integral  inversions,  that  minimized  the 
discrepancies  between  the  input  curves  with  back  calculated  results,  down  to  the  level  of  the  scatter 
of  the  reduced  curves  (Rheological  Transforms  Software  RTS:  Time-Temperature  Instruments,  P. 
O.  Box  40156,  Pittsburgh,  PA  15201). 

The  real  G1  and  imaginary  G"  components  of  the  complex  dynamic  shear  modulus  G*  =  G1 
+  iG//  for  the  nascent  molecular  network  (Epon  1001F/0.55HA/0.45DDS)  are  presented  in  Fig.  15 
along  with  a  linear  plot  of  the  loss  tangent.  Tan  5  *  Gl/G/  =  J#/J/,.  Except  in  the  region  of  glassy 
response  at  high  reduced  frequencies,  wr  =  waT,  the  tenuous  gel  is  very  lossy  with  the  loss  tangent  with 
values  which  exceed  two.  It  is  clear  that  exponential  behavior  is  seen  between  reduced  frequencies 
of  10‘8  and  10' 13  where  the  slopes  of  log  G/  and  log  G#  are  both  about  0.85.  If  true  exponential 
behavior  were  displayed  the  loss  tangent  would  be  constant  in  the  same  frequency  range.  The  low 
frequency  peak  is  rather  broad.  The  absence  of  constancy  may  be  due  to  the  fact  that  this  sample 
is  very  slightly  beyond  the  gel  point. 

To  illustrate  how  the  inclusion  of  the  non-additive  effects  of  the  viscosity  and  the  glassy 
compliance  in  the  modulus  destroys  the  systematic  and  readily  interpretable  changes  seen  in  the 
compliance  curves,  most  of  the  curves  seen  in  Fig.  8  have  been  transformed  to  G^w)  and  G^(w). 
The  results  are  shown  in  Figs.  16  and  17  respectively.  The  low  frequency  limits  of  Log  G\u)  for 
0.50, 0.70  and  1 .0  DDS  are  the  equilibrium  moduli  Ge  of  the  increasingly  tighter  molecular  networks. 
The  limiting  low  frequency  slopes  of  log  Gf(u>)  and  log  Gf(w)  for  the  nascent  molecular  network 
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sample  0.45  DDS  are  both  0.85.  The  low  frequency  slopes  of  log  G^w)  of  the  viscoelastic  liquid 
samples  0.25,  0.35  and  0.40  DDS  are  2.0,  1.4,  and  1.7.  The  0.25  DDS  sample  has  reached  the  limiting 
value  of  2.0  but  the  0.35  and  0.40  DDS  samples  have  not.  Presumably  when  they  do  they  will  be  in 
inverse  order  on  the  log  ui  scale  according  to  their  concentration  of  DDS;  i.e.  the  0.25  DDS  curve  will 
be  found  at  the  highest  frequencies  with  the  0.35  and  0.40  DDS  curves  at  correspondingly  lower 
frequencies.  At  the  lowest  frequencies  in  Fig.  17  it  can  be  seen  that  log  G#(w)  curves  foi  all  of  the 
materials  with  the  exception  of  the  0.45  DDS  sample  have  reached  the  limiting  slope  of  one.  For  the 
viscoelastic  liquids  this  represents  complete  domination  by  viscous  flow  since  in  this  region  Gf(w)  = 
wij,  where  tj  is  the  steady  state  shear  viscosity  coefficient.  For  the  viscoelastic  solids  with  microscopic 
molecules  networks  the  real  component  of  the  dynamic  viscosity  has  reached  a  constant  low 
frequency  limit.  With  both  sets  of  modulus  curves  we  find  very  little  that  is  apparent  or  interpretable 
in  terms  of  the  changing  molecular  structure  on  the  molecular  network  is  developed.  The  GA(w)  are 
indeed  the  most  unenlightening.  Finally  the  corresponding  relaxation  spectra  for  the  different  stages 
in  network  development  are  shown  in  Fig.  18.  We  find  them  to  be  an  unintelligible  morass  of 
crossing  curves. 
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TABLE  II 


Viscosities  of  Non-gelled  Epons 


log  rj  = 

log  A 

DDS  Fraction  log  A 

C 

0.25 

-5.14 

1850 

0.35 

-3.419 

1850 

0.40 

-2.670 

1850 

C/CT-T J  ] 

T. 

log  Tj(T0  J) 

•c 

poise 

15.8 

11.49 

10.0 

10.93 

10.0 

10.74 

LEGENDS 


Fig.  1.  Logarithmic  presentation  of  the  recoverable  shear  compliance  Jr(t)  of  Epon  100 IF  as  a 
function  of  the  logarithm  of  time  t  at  9  temperatures  as  indicated. 

Fig.  2.  Logarithmic  plot  of  the  steady  state  recoverable  shear  compliance  Je°  of  the  neat  resin  EPON 
100 IF  as  a  function  of  temperature. 

Fig.  3.  Attempt  to  reduce  the  recoverable  compliance  curves  to  a  common  time-scale  of  response  by 
using  the  viscosity  ratio  »KT0)/»?(T)  as  the  temperature  shift  factor.  The  chosen  references 
temperature  T0  =  30.2°C. 

Fig.  4.  Superposition  failure  for  the  Epon  1001F  data  from  three  temperatures  shown  by  plotting  the 
reduced  recoverable  compliance,  Jr  ,  as  a  function  of  the  reduced  time-scale. 

Fig.  5.  Logarithm  of  the  recoverable  compliance  Jr(t)  of  the  EPON  100 IF  epoxy  resin  fully  reacted 
with  0.25  equivalent  ratio  of  DDS  and  0.75  of  MA  plotted  against  the  logarithm  of  the 
reduced  time,  t/ar  Recoverable  compliance  curves  from  five  temperatures  were  reduced 
to  a  reference  temperature  of  64.0°C.  The  dashed  line  represents  the  creep  compliance  curve, 
J(t). 

Fig.  6.  Reduced  log  Jr(t)  vs  Log  t/aT  plot  for  epoxy  resin  Epon  100 IF  reacted  with  35%  crosslinking 
amine,  DDS.  Five  curves  determined  at  the  indicated  temperatures  were  used  to  produce  the 
curve  reduced  to  65°C  as  the  reference  temperature  T0. 
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Fig.  7.  Reduced  compliance  Jp(t)  for  epoxy  cured  with  45%  DDS  presented  in  similar  fashion  to  the 
J(T)  curves  in  Figs.  3  and  4.  Six  curves  were  shifted  to  72.50C.  This  material  was  exceedingly 
close  to  the  incipient  point  of  gelation. 

Fig.  8.  Comparison  of  all  of  the  time  dependent  reduced  compliance  [Jp(t)  -  Jg]  dashed  line  curves 
as  logarithmic  functions  of  log  t/aT.  Comparison  temperatures  are  indicated  above  and  in 
Table  I. 


Fig.  9.  Double  logarithmic  plot  of  the  retardation  spectrum  L(r)  as  a  function  of  the  reduced 
retardation  time  obtained  from  the  data  of  Fig.  3  for  the  epoxy  cured  with  0.25%  DDS. 

Fig.  lO.Double  logarithmic  plot  of  L(r)  shown  as  a  function  of  r  obtained  from  the  data  of  Fig.  4  for 
an  epoxy  cured  with  0.35%  DDS. 

Fig.  11. Double  logarithmic  comparison  plot  of  all  of  the  L(r)s  obtained  on  the  fully  reacted  Epon 
1001F/DDS/MA  mixtures.  The  dashed  lines  indicate  the  absence  of  a  macroscopic  network 
(gel).  The  solid  lines  indicate  that  a  molecular  network  is  present. 

Fig.  12.VFTH  plot  of  the  time  scale  shift  factors  of  aT  obtained  from  the  reduction  process  to  obtain 
the  compliance  curves  of  Fig.  6.  See  text  for  details. 
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Fig.  13.Logarithmic  plot  of  the  shear  creep  compliance  measured  during  the  curing  of  Epon  100 IF 
with  a  stoichiometric  amount  of  DDS  at  142. 8°C.  All  measurements  were  made  after  gelation 
occurred.  Curing  time  Tc  when  creep  measurements  were  made  are  indicated.  Long-dashed 
extrapolation  lines  calculated  as  described  in  text.  The  short-dashed  lines  are  empirically 
shifted  to  make  reasonable  extrapolations  during  the  latter  stages  of  cure. 

Fig.  14.  Selected  curves  from  Fig.  13  reduced  to  the  curve  with  tc  =  22.4  hrs.  of  cure,  after  the 
subtractions  of  the  glassy  compliance  Jg. 

Fig.  15.The  storage  G1  and  loss  G 1  moduli  of  the  incipient  gel  (Epon  1001/0.55  MA/0.45  DDS) 
plotted  logarithmically  as  a  function  of  the  logarithm  of  the  reduced  frequency  =  waT 
(radius  sec'1).  The  loss  tangent  is  also  shown.  The  temperature  of  reduction  is  12.5°C. 

Fig.  16.The  storage  modulus  G7  plotted  logarithmically  as  a  function  of  the  reduced  frequency  for 
the  different  stages  of  network  development  of  the  Epon  1001F/DDS  system.  A  solid  line 
indicates  a  viscoelastic  solid.  Long-dashed  lines  indicate  viscoelastic  liquids.  Short-dashed 
lines  are  extrapolations  beyond  the  reduced  measurement  range.  The  reference  temperatures 
are  the  same  as  those  indicated  in  Fig.  8. 

Fig.  17.The  loss  modulus  G 1  for  the  various  stages  in  network  development  plotted  according  to  the 
description  of  Fig.  16. 

Fig.  18.The  relaxation  spectra  H(r)  for  the  different  stages  of  network  development  presented  as 


described  in  Fig.  16. 


Log  Jr(t) 


Figure  3 


Figure  4 
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Figure  5 


Figure  6 


LogJ^(t)  Pa'1 


Figure  12 
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Figure  13 


Figure  14 


Figure  15 
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XII.  EPOXY  RESINS  (DGEBA),  the  curing  and  physical  aging  process 

By 

D.  J.  Plazek*  and  Z.  N.  Frund,  Jr.** 

ABSTRACT 

Differential  Scanning  Calorimetry  (DSC)  and  Infrared  Spectroscopy  (IR) 

were  employed  to  monitor  the  degree  of  cure  of  partially-cured  epoxy  resin 

(Epon  828/MDA)  samples.  The  extent  of  cure,  as  determined  by  residual  heats 

of  reaction  concurred  with  that  determined  by  monitoring  the  infrared 

radiation  absorbance  of  the  epoxide  group  near  916  cm  .  The  fictive 

temperature  T^  was  found  to  increase  with  the  degree  of  cure,  increasing 

rapidly  during  cure  until  reaching  a  value  near  the  cure  temperature  T^  of 

130°C  (approximately  80  percent  cure)  where  the  material  vitrified.  The 

greatly  reduced  reaction  rate  during  the  final  20  percent  of  cure  was  not  only 

a  consequence  of  vitrification  but,  as  revealed  by  infrared  spectroscopy,  the 

result  of  the  depletion  in  the  number  of  reactive  epoxide  groups.  The 

endothermic  peak  areas  and  peak  temperatures  evident  during  the  DSC  scans  were 

used  as  a  measure  of  the  extent  of  "physical  aging"  which  took  place  during 

the  cure  of  this  resin,  and  after,  fully-cured  samples  were  aged  37  C  below 

their  ultimate  glass  temperature  for  various  periods  of  time.  The  rate  of 

physical  aging  slowed  down  as  the  temperature  increment  (T^  -T^)  increased. 

While  an  endothermic  peak  was  evident  after  only  one  hour  of  cure  (T^  = 

138. 3°C),  such  as  a  peak  did  not  appear  until  fully-cured  samples  were  aged 

for  16  hours  or  more.  Enthalpy  data  revealed  that  for  partially-cured 

material,  the  fictive  temperature  Tf  ,  reflecting  physical  aging,  increased 

with  curing  time.  In  contrast,  the  Tf  for  fully-cured  samples  decreased 

with  sub-T  aging  time.  The  characteristic  jump  in  the  heat  capacityACp 

which  occurred  at  the  T,  decreased  as  curing  progressed.  This  decrease  appears 

*  |  ^  _ 
to  be  dependent  upon  the  rotational  and  vibrational  degrees  of  freedom  of 

the  glass.  Finally,  a  graphical  method  of  determining  the  fictive  temperature 

T,  of  partially  and  fully-cured  epoxy  material  from  measured  endothermic 
*  >  ® 

peak  areas  was  developed  . 

*Unlversity  of  Pittsburgh,  Department  of  Material  Science  and  Engineering, 
Pittsburgh,  PA 

**Mine  Safety  Appliances  Company,  Pittsburgh,  PA 
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INTRODUCTION 


Epoxy  resins,  traditionally  used  only  as  adhesives  and  coatings  have 

become  the  most  common  matrix  for  aerospace  composite  materials.  The  degree 

of  cure  of  an  epoxy  is  influential  in  determining  its  physical  properties, 

i.e.,  dimensional  and  thermal  stability,  compliance,  glass  temperature,  etc. 

Of  particular  concern  are  the  consequences  of  physical  aging  (densif ication  at 

a  temperature  below  the  glass  temperature  T  )  (1,2),  one  of  which  is  the 

§ 

increase  in  the  material's  brittleness. 

Previous  studies  (3,4)  have  been  made  on  the  physical  properties  of 
epoxy  resins  during  and  after  curing.  During  curing  the  increasing  density 
was  monitored  in  dilatometer.  The  increasing  density  was  caused  by  the 
substitution  of  shorter  primary  covalent  bonds  for  secondary,  dispersion  or 
van  der  Uaals  "bonds"  as  the  chemical  reaction  proceeded.  The  higher 
densities  reflect  greater  molecular  crowding;  both  of  which  are  manifested  by 
the  higher  glass  temperatures  T 


It  was  observed  that  when  the  T  of  the  system  became  greater  than  the 

8 

curing  temperature  Tc>  not  only  did  the  reaction  deccelerate  but  endothermic 
peaks  in  the  heat  capacity  -  temperature  curves  and  maxima  in  the  thermal 
expansion  coefficients  c<  appeared  and  increased  with  time.  It  became  clear 
that  these  peaks  were  simply  manifestations  of  what  has  become  known  as 
"physical  aging". 

When  an  amorphous  material  is  held  at  a  temperature  below  its  T  ,  the 

8 

density  increases  spontaneously  with  the  passage  of  time  as  an  equilibrium 
liquid  structure  is  approached.  In  addition,  many  of  its  kinetic  responses 
become  sluggish.  Upon  heating,  such  an  aged  glass  retains  the  glassy  values 
for  the  Cp  andc<  to  temperatures  above  T^.  As  the  material's  temperature 


sp 


increases  further,  acceleration  in  the  enthalpy  H  and  specific  volume  v 
curves  are  observed.  As  higher  molecular  mobilities  are  achieved  during 
heating  the  material  can  reach  an  equilibrium  structure.  This  increase  in 
mobility  is  autocatalytic  since  the  increasing  free  volume  which  occurs  with 
thermal  expansion  is  the  rate  determining  parameter. 

Once  the  increases  past  Tc  the  material  is  in  the  physical  aging 
domain.  Network  development  and  physical  aging  are  concurrent.  The 
endothermic  peaks  found  at  this  stage  of  curing  were  so  large  that  it  was 
thought  that  densif ication  was  more  efficient  during  reaction  than  after  full 
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curing.  It  was  surmised  that  the  mobility  of  the  unattached  molecules  was 

possibly  greater  than  that  of  the  molecular  segments  incorporated  into  the 

network.  A  study  of  the  aging  which  occurred  during  the  curing  reaction  and 

after  full  cure  would  have  to  be  performed  to  study  this  hypothesis.  Such  a 

study  was  carried  out  and  is  reported  here. 

Differential  Scanning  Calorimetry  (DSC)  and  Infrared  Spectroscopy  (1R) 

were  used  to  monitor  the  heats  of  reaction,  and  infrared  absorbance, 

respectively,  of  a  diglycidyl  ether  of  bisphenol  A  (DGEBA)  epoxy  resin  cured 

for  different  periods  of  time.  Differential  Scanning  Calorimetry  was  also 

used  to  measure  the  extent  of  physical  aging  which  took  place  during  the  cure 

of  an  epoxy  resin  at  a  temperature  below  T  (oO)  the  glass  temperature  of  the 

S 

fully-cured  resin.  Fully-cured  epoxy  samples  were  also  aged  at  the  chosen 
curing  temperature. 


BACKGROUND 


The  Curing  Process 

Curing  of  an  epoxy  resin  involves  the  development  of  a  three-dimensional 
molecular  network  from  the  reactants  via  branching  during  polymerization  with 
no  volatile  by-products.  The  combination  of  an  epoxide  and  the  primary  amine 
leads  to  two  principal  reactions.  The  first  is  the  expected  attack  of  the 
epoxide  group  by  the  primary  amines.  When  the  remaining  secondary  hydrogens 
combine  with  the  second  epoxy  molecule  a  branch  point  is  formed.  As  a  result 
of  branching,  thermosetting  resins  become  set,  i.e.,  infusible  and  insoluble. 
Whether  branching  or  linear  growth  occurs  faster  depends  on  the  relative 
reaction  rate  of  the  epoxide  with  the  primary  or  secondary  amine  (5). 
Eventually  the  molecular  chains  become  linked  together  into  networks  of 
effectively  infinite  molecular  weight.  The  volume  contraction  which  occurs 
during  cure  is  a  result  of  the  exchange  of  van  der  Waals  "bonds”  for  shorter 
covalent  bonds.  Because  of  increased  molecular  crowding,  the  glass 
temperature  is  elevated. 
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A  1  ?H 

H2Nwww  NH2  +  H2C“CH  vaw — >H2NwwN  —  CHj— CHwww 

Primary  Amine  Epoxide  Primory,  Secondory  Amines  ond 

Pendant  Hydroxyl 

H  OH  O  OH 

HNvkwN-  CH-CH  WWW+  HC-CH  wwv-oHNvww  N-CH  -  CH  vvwvw 

*  *  1  *  |  2 

Rrimory  and  Secondary  Epoxide  CH2  Tertiary  Amines 

HC-OH 

Figure  1.  Chemical  Reaction  of  Primary  and  Secondary  Amines  with  Epoxide 

Group 

•The  moment  of  a  sudden  and  irreversible  transformation  from  a 
viscoelastic  liquid  to  a  solid  gel  is  known  as  the  gel  point.  At  this  point  a 
covalently  bonded  molecular  network  extends  throughout  and  the  viscosity 
becomes  infinite.  Gelation  occurs  at  a  specific  point  in  the  course  of  the 
chemical  reaction  and  depends  on  the  functionality,  reactivity,  and 
stoichiometry  of  the  reactants  (6).  Gelation  typically  occurs  between  55  and 
80  percent  conversion  (degree  of  cure  p  =  0.55  -  0.80). 

Distinct  from  gelation,  vitrification  may  occur  at  any  stage  of  the 
reaction.  This  transformation  from  a  viscous  liquid  or  elastic  gel  to  a  glass 
appears  when  the  glass  transition  temperature  of  the  material  reaches  the  cure 
temperature.  Further  cure  in  the  glassy  state  is  slow  and  usually  diffusion 
controlled. 

Vitrification  is  a  reversible  process.  Curing  can  be  accelerated  by 
heating  the  partially-cured  thermoset  above  its  glass  temperature  T^. 
Therefore,  to  assure  full  cure,  thermosets  should  be  cured  near  or  above  their 
ultimate  glass  temperature. 

The  Glass  Temperature  T^  Versus  Flctive  Temperature  T^ 

The  glass  temperature  T^  is  the  temperature  below  which  a  material 
exhibits  properties  typical  of  a  glass.  It  is  usually  defined  in  terms  of  the 
material's  specific  volume  or  enthalpy  as  a  function  of  temperature,  obtained 
at  a  specified  cooling  rate.  At  temperatures  above  the  glass  transition  the 
molecular  mobility  is  great  enough  for  the  liquid  to  rearrange  rapidly  to  a 
structure  characteristic  for  its  temperature.  Upon  cooling,  the  increase  in 
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molecular  crowding  decreases  the  mobility  to  the  point  where  rearrangements 
necessary  to  achieve  an  equilibrium  liquid  configuration  and  density  cannot 
keep  up  with  the  rate  of  cooling.  The  point  of  departure  from  this 
equilibrium  state  or  break  (abrupt  decrease  in  slope)  in  the  specific 
volume-temperature  cooling  curve  depends  on  the  rate  of  cooling.  The 
departure  or  break  comes  sooner  -  at  a  higher  temperature  with  faster  cooling. 
At  slower  cooling  rates  the  sluggish  molecular  rearrangements  can  keep  up  with 
cooling  until  a  lower  temperature  is  reached.  This  break  point  is  identified 
as  a  glass  transition  temperature  and  is  a  material  characterizing  parameter. 

According  to  Toole  (7),  the  fictive  temperature  is  a  hypothetical 
temperature  at  which  the  structure  of  the  glass  would  be  at  equilibrium.  It 
can  be  determined  from  DSC  heating  curves.  The  response  of  any  glass  upon 
heati'ng  is  a  function  of  its  thermal  history  (i.e.,  time  of  aging  at  a 
constant  temperature,  rate  of  cooling,  etc.)  and  the  rate  of  heating.  The 
fictive  temperature  is  a  specimen  property  rather  than  a  characteristic 
material  property.  Although  Toole's  hypothesis  may  not  be  correct,  the 
operational  definition  is  useful  in  characterizing  glassy  specimens.  The 
fictive  temperature  can  be  determined  by  the  intersection  point  of  the 
measured  specific  volume~(or  enthalpy)  temperature  line  with  the  equilibrium 
volum.e-(or  enthalpy)  temperature  line  of  the  specimen. 

The  glass  transition  temperature  must  be  determined  from  cooling,  not 

heating  curves.  Unfortunately,  because  of  unavoidable  supercooling  of 

calibration  standards,  accurate  T  values  cannot  be  measured  with  current 

8 

differential  scanning  calorimeters  since  these  instruments  cannot  be  corrected 
for  temperature  lags.  Even  if  these  lags  could  be  corrected  for, 
partially-cured  thermosets  cannot  be  taken  to  temperatures  which  are  much 
higher  than  the  curing  temperature  (to  prepare  for  cooling)  without  incurring 
further  chemical  reaction.  This  additional  reaction  changes  the  value  being 
sought.  Consequently,  only  the  fictive  temperatures  of  curing  systems  can  be 
determined. 

In  spite  of  the  fact  that  glass  transition  temperatures  can  only  be 

determined  from  cooling  experiments,  most  reported  "glass  temperatures"  are 

obtained  from  heating  curves  or  from  measurements  of  loss  moduli  or  loss 

tangents.  The  latter  qualitatively  correlate  with  T  's  but  are  always  higher 
o  8 

by  15  or  20  C.  The  chosen  frequencies  determine  the  position  of  the 

mechanical  property  maxima.  In  a  material  with  minimal  aging  the  T^  will  not 

be  more  than  a  few  degrees  below  the  T  at  the  same  rates  of  heating  and 
cool ing. 
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MATERIALS 


Epoxy  resir.s  are  characterized  by  a  chemical  structure  containing  the 
highly  reactive  epoxide  group  or  oxirane  ring. 


Hx  'H 


Figure  2.  Chemical  Structure  of  Epoxide  Group 


'  In  the  present  study,  the  Epon  828  epoxy  resin  was  cured  with  methylene 
dianiline  (MDA).  Epon  828,  principally  a  diglycidyl  ether  of  bisphenol  A,  is 
difunctional.  The  tetrafunctional  methylene  dianiline  is  produced  by 


Uniroyal.  The  chemical  structure  and  some  characteristics  of  the  Shell  Epon 
828  and  Uniroyal  Tonox  are  shown  below. 

Shell  Epon  828  Epoxy  (DGEBA)  Resin  Structure  and  Characteristics 

CH, 

'CH-  CHjOO-  C-O-  0-  CH  —  CH-'cH 
CH, 


h2c 


2 


Viscosity  (25  C) :  110-150  poise 

Specific  Gravity  (20°C):  1.10 


Epoxide:  0.52-0.54  eq./llO  gms 
M  :  380  gms/mol 

n  T  :  -14°C 

Epoxy  Equivalent®Weight  (number  of  grams 
of  resin  containing  one  chemical  equivalent 
of  an  epoxy  group):  190 


Functionality  (M^/epoxy  equivalent  weight): 


2.0 


Uniroyal  Tonox  Methylene  Dianiline  (MDA) 


h2n-0-  ch2-0-nh2 

Specific  Gravity  (20°C):  1.16 

Molecular  Weight:  198.3  cms/mol 
Melting  Point:  90  ‘3  C 
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EXPERIMENTAL  PROCEDURES 

Stoichiometric  ratios  of  the  Epon  828  and  MDA  were  mixed.  The  mixture 
was  heated  in  a  Pyrex  beaker  to  95°C  and  stirred  by  hand  until  a  clear 
homogeneous  solution  resulted.  This  beaker,  with  the  mixture,  was  then 
partially  immersed  in  ice  water  for  30  minutes  to  inhibit  further  reaction. 
Heats  of  reaction  and  thermal  phenomena  near  the  glass  temperature  T^  were 
measured  with  a  Perkin-Elmer  System  4  Differential  Scanning  Calorimeter.  An 
Analect  AQS-20  Fourier  Transform  Infrared  (FTIR)  Spectrophotometer  was 
employed  to  monitor  the  changes  which  occurred  within  the  absorbance  band 
intensities  of  partially-cured  material. 


Heats  of  Reaction  and  Thermal  Phenomena  During  Cure 

.  Endothermic  peaks  reflecting  the  state  of  densif ication  of  aged  samples 

were  determined  for  samples  at  different  stages  of  cure.  This  measurement 

erased  the  effects  cf  aging  and  the  fictive  temperature  of  the  samples  were 

then  determined  immediately  after  cooling. 

For  the  DSC  investigation,  9-14  milligrams  of  the  uncured  mixture  were 

measured  into  standard  aluminum  pans  and  covered  with  lids.  Samples  were  then 

placed  in  0.625  inch  diameter  by  0.375  inch  deep  cavities  within  a 

thermostatted  metal  block.  The  metal  block  was  preheated  and  maintained  at 
o 

130+_1  C.  Samples  were  periodically  removed  from  the  metal  block  and  placed  on 
a  metal  surface  at  23°C.  Within  30  minutes  the  endothermic  peak  areas  and 
temperatures,  and  fictive  temperatures  were  determined  at  a  heating  rate  of 
20  C/minute.  The  temperature  and  heat  flow  calibration  of  the  DSC  was 
performed  with  indium  and  synthetic  sapphire  standards  at  different  heating 


To  determine  the  degree  of  cure  samples  were  removed  from  the  heated 
metal  block  and  placed  in  the  DSC  at  0°C  where  negligible  reaction  took  place. 
They  were  then  heated  to  169°C  at  a  rate  of  100°C/minute.  Samples  were 
maintained  at  169°C  for  as  long  as  one  hour.  The  total  exotherm  during  this 
period  was  used  as  a  measure  of  the  samples'  degree  of  cure.  At  169°C  all 
samples  were  found  to  achieve  their  ultimate  fictive  temperature  of  167°C 
within  one  hour  or  less.  A  heating  rate  of  100°C/minute  was  used  to  minimize 
any  unrecorded  reaction. 
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Infrared  Spectroscopy  During  Cure 

The  absorption  of  infrared  radiation  due  to  the  molecular  vibration  of 
epoxide  groups  at  916  cm  *  was  measured  as  a  function  of  cure  time.  A 
microthin-layer  of  the  uncured  mixture  was  evenly  spread  between  two  potassium 
bromide  (KBr)  windows.  Prior  to  this  however,  the  windows  were  polished  and 
their  spectra  obtained  and  stored  in  the  FTIR  spectrophotometer  as  a 
background  reference.  The  KBr  windows  with  sample  were  placed  between  two 
metal  plates  at  130+1°C  for  discrete  periods  of  time.  They  were  then  placed 
in  the  spectrophotometer  and  a  spectrum  (4000-400  cm  * )  from  32  scans 
obtained. 


Thermal  Phenomena  During  Physical  Aging  With  No  Reaction 
•To  monitor  the  physical  aging  process  without  a  simultaneous  chemical 
reaction,  fully-cured  DSC  samples  in  which  the  previous  effects  of  aging  had 
been  erased  were  placed  back  into  the  heated  metal  block  at  130j+1°C. 
Fully-cured  samples  were  used  so  that  the  aging  process  was  not  masked  or 
effected  by  the  curing  reaction.  Periodically,  samples  were  removed  from  the 
metal  block  and  quickly  cooled.  The  endothermic  peak  temperature  and  area, 
and  fictive  temperature  of  these  samples  were  then  determined. 

RESULTS  AND  DISCUSSION 

Thermal  Phenomena  and  Heat  of  Reaction  During  Cure 
-  Table  I  is  a  listing  of  the  thermal  phenomena  and  degree  of  cure  values 
determined  from  DSC  studies  and  infrared  spectroscopy  for  partially-cured 
samples.  After  curing  for  one  hour  or  more  the  fictive  temperatures  reached 
values  greater  than  the  cure  temperature  Tc  of  130°C  and  a  dense  glass  was 
formed.  This  was  a  simple  consequence  of  "physical  aging"  during  cure  and  is 
manifested  by  a  large  endothermic  peak  upon  first  DSC  heatings.  See  Figure  3. 
Sickfeld  and  Mielke  (8)  reported  the  presence  of  an  endothermic  peak  when 


Figure  3  revealed  that  the  thermograms  had  similar  slopes  and  that  peak 
areas  and  temperatures  increased  with  cure  time.  Such  samples  were  quickly 
cooled  at  80°C/minute  from  just  beyond  their  endothermic  peak.  This  brief 
excursion  above  the  glass  temperature  erased  the  effects  of  physical  aging  yet 
did  not  cause  appreciable  additional  cure.  A  second  heating  which  immediately 
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followed  was  necessary  to  obtain  the  "normally  shaped"  curves  presented  in 
Figure  4. 

Two  sets  of  fictive  temperatures  are  listed  as  a  function  of  curing  time 

in  Table  I.  The  lower  values,  T,  ,  were  determined  from  enthalpy  values,  H, 

1 9  a 

obtained  during  first  DSC  heating  scans.  They  reflect  the  amount  of  physical 
aging  which  occurred  during  the  curing  reaction,  after  the  material's  T 

8 

became  greater  than  the  curing  temperature.  In  Figure  5,  the  AH  values  for  a 

sample  cured  for  nine  hours  at  130_+1°C  are  plotted  as  a  function  of 

temperature.  The  Ah  values  were  determined  from  the  equation 

AH=H(T)  -  [H(150°C)  +  C  (T  -  150°C)],  where  the  heat  capacity  C  was 

Pg  Pg 

determined  from  the  slope  of  the  enthalpy  H  versus  temperature  curve  in  the 

glassy  region.  A  T,  of  151. 0°C  was  found  by  extrapolating  back  from  the 
*  > a  0 

linear  portion  of  the  curve  (above  180  C  in  this  case)  to  the  portion  of  the 

curve  having  a  slope  close  to  zero.  Plots  of  AH  versus  T  showed  such  slight 

changes  in  slope  that  extrapolations  were  not  definitive.  The  correspondingly 

higher  T,  values  were  obtained  directly  from  the  second  heating  curves  after 
‘  t ,  g  ~ 

aging  had  been  erased  by  the  initial  heating  to  a  temperature  above  T  . 

g 

As  illustrated  in  Figure  4,  the  fictive  temperature  Tr  was  taken  as 

f » g 

the  midpoint  from  the  extension  of  the  pre-  and  post-transition  baselines  of 

the  DSC  thermogram.  The  T^  is  the  point  in  the  thermogram  where  the  heat 

capacity  changes  the  most  rapidly.  All  thermograms  had  similar  slopes  and 

their  T,  's  increased  with  cure  time.  It  can  be  seen  in  Figure  4  that  when 
*  8 

there  is  little  or  no  physical  aging,  there  is  no  maximum  in  the  heat  capacity 

curve.  A  sample's  T^  is  expected  to  be  close  to  its  glass  temperature 

(normally  obtained  during  cooling)  if  the  heating  follows  immediately  after 

cooling  from  above  T  (3). 

g 

The  total  peak  areas  values  listed  in  Tables  I  and  II  were  determined  by 
superimposing  the  second  DSC  heating  scan  of  a  sample  over  its  first.  The 
peak  area  bounded  by  these  curves  was  then  measured.  This  is  illustrated  in 
Figure  3  for  a  sample  cured  for  336  hours  at  130°C. 


A  sample's  degree  of  cure  p  was  calculated  from  the  equation 

p  =  (Ah  -  Ah  )/AH  where  AH  is  the  total  exotherm  for  the  initially 
o  c  o  o 

uncured  sample  and  AH^  is  the  exotherm  of  a  partially-cured  sample.  The  Epon 
828/MDA  was  heated  to  169°C  immediately  after  mixing  (with  no  induced  cure). 
Its  exotherm  yielded  a  AHq  equal  to  98.7  calories/gram.  Figure  6  presents 
the  isothermal  DSC  thermograms  for  an  uncured  sample  and  those  which  had  been 
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partially-cured  for  0.166  and  0.250  hours  at  130+^1°C.  It  illustrates  how  the 
exotherm  decreased  when  a  sample  was  cured  for  increasing  periods  of  time. 
After  approximately  40  minutes  at  169°C,  all  samples  achieved  the  ultimate 
T,  of  167°C. 

f’8  o  o 

The  fact  that  perceptible  curing  proceeds  near  130  C,  a  temperature  37  C 

below  the  material's  ultimate  T c  where  molecular  crowding  is  so  great  that 

the  equilibrium  density  cannot  be  reached  in  less  than  a  millenium  is 

suprising.  Aherne,  Enns,  Doyle  and  Gillham  (9)  have  reported  similar  findings 

for  epoxy  when  T  =  T  -  26°C. 
v  1  eg 

Infrared  Absorbance  During  Cure 

The  infrared  spectra  of  this  epoxy  resin/MDA  system  exhibited  an 
absorption  band  due  to  the  vibration  of  the  epoxide  groups  near  916  cm  See 

Figure  7.  This  vibration  phenomenon  was  due  to  the  bending  of  the  -C-0-C- 
triad.  According  to  Antoon,  Starkey  and  Koenig  (10)  and  others  (11),  an 
epoxy's  extent  of  reaction  can  be  monitored  by  measuring  the  intensity  of  the 
916  cm  £  absorption  band,  Ag1(.  cm_£ *  Figure  7  illustrates  how  the  intensity 
of  the  916  cm  1  band  decreases  with  increasing  cure  time,  indicating  a 
reduction  in  the  number  of  epoxide  groups.  The  1-Ag^  cra-1  values  for 
partially-cured  samples  are  listed  in  Table  I. 

Bands  initially  at  3340  and  3473  cm-1  arose  from  the  stretching  of  the 
N-H  bonds  within  the  MDA.  As  curing  progressed,  the  amine  hydrogens  were 
depleted  and  a  broad  0-H  stretch  between  3200  and  3600  cm  became  evident. 
This  was  due  to  the  hydrogens  of  the  MDA  bonding  to  the  oxygens  from  the 
reacted  epoxide  groups.  As  curing  progressed,  the  quantity  of  -OH  groups 
increased.  The  0-H  band  was  very  broad  because  of  hydrogen  bonding  (12). 


Kinetics  of  Curing 

Three  different  measured  quantities  -  the  degree  of  cure  obtained  from 

hcatr.  of  reaction,  the  fictive  temperature  T,  and  the  infrared  absorbance  of 
-1  * 

the  916  cm  band,  listed  in  Table  I  were  plotted  as  a  function  of  the 

logarithmic  time  of  cure.  See  Figure  8.  Within  experimental  uncertainty,  the 

data  points  for  all  three  quantities  fell  on  the  seme  sigmoidal  curve.  It  was 

shown  by  Choy  and  Plazek  (3)  that  the  T,  monitored  the  extent  of  cure  which 

£  >  8 

had  been  calculated  from  the  volume  shrinkage  during  the  cure. 

When  plotted  against  log  t  ,  the  fictive  temperature  T,  of  the  Epon 
e  °  cure  r  f,g 
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828/MDA  rapidly  increased  until  reaching  a  value  near  130  C  where  it  began  to 

vitrify.  Subsequently,  an  abrupt  leveling  out  was  observed  as  the  Tf 

*  *  8 

exceeded  the  cure  temperature  Tc«  Ultimately,  the  material  achieved  a  T, 

o  f  ® 

value  near  167  C. 

The  degree  of  cure  (DSC)  and  l-A^^  -1  values  rapidly  increased  until 

reaching  approximately  80  percent  of  total  cure.  The  greatly  reduced  rate  of 
reaction  during  the  final  20  percent  of  cure  was  not  only  a  consequence  of 
vitrification  but,  as  revealed  by  infrared  spectroscopy,  the  result  of  a 
depletion  in  the  number  of  reactive  epoxide  groups.  The  absorbance  reached  an 
asymptotic  limit  of  0.06.  This  may  reflect  isolated  epoxide  groups  which 
cannot  reach  amine  hydrogens  and  react  in  any  practical  time  at  any 
temperature,  c.r  reflectance  of  infrared  radiation  at  the  surface  of  the  KBr 
window. 


Physical  Aging  Phenomenon  Without  Simultaneous  Reaction 
The  extent  of  physical  aging  at  130°C  was  determined  for  fully-cured 
Epon  828/MDA  samples  aged  for  as  long  as  648  hours.  As  with  the  curing 
process,  high-dersity  glasses  were  formed.  Such  densified  glasses  responded 
sluggishly  to  temperature  changes  near  the  fictive  temperature.  Figure  9 
illustrates  that  upon  first  heatings,  the  DSC  thermograms  had  similar  slopes 
and  that  the  endothermic  peak  temperatures  and  areas  increased  with  the  time 
of  aging.  Figures  10  and  11  revealed  that  both  the  peak  temperatures  and 
areas  were  linear  functions  of  the  log  tag^ng*  Kong  (13)  previously  found 
that'a  linear  relationship  existed  between  the  endothermic  peak  areas  and 
aging  time. 

Although  an  endothermic  peak  was  evident  after  only  one  hour  of  cure 

(T^  g=138.3°C),  such  a  peak  did  not  appear  until  fully-cured  samples  were  aged 

for  16  hours  or  more.  Physical  aging  takes  place  faster  when  a  sample  is  aged 

at  a  temperature  nearer  to,  yet  below,  its  fictive  temperature.  Therefore,  it 

is  not  surprising  that  the  effects  of  physical  aging  were  evident  earlier  and 

in  greater  magnitude  in  partially-cured  samples  having  a  T,  closer  to  the 

* » 8 

aging  temperature.  Kong  (13)  demonstrated  that  the  kinetics  of  aging  within 

epoxy  resins  slowed  down  as  the  increment  (T  -  T  )  increased. 

8  a 

Sickfeld  and  Heinze  (14)  observed  endothermic  peaks  only  after  an  epoxy 
resin's  T  exceeded  its  annealing  temperature  by  at  least  10°K.  As  in  this 

O 

investigation,  the  temperature  of  the  peaks  increased  with  aging  time. 
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Sickfeld  surmised  that  the  endothermic  peaks  were  due  to  a  higher  degree  of 
order  in  the  glassy  state. 

Kreibich  (15)  evaluated  the  effects  of  aging  in  hexahydrophthalic  acid 

glycidyl  ester/hexahydrophthalic  acid  anhydride  and  found  the  peak  areas  to 

also  increase  with  time.  She  concluded  that  by  aging  the  material  slightly 

below  its  the  polymer  chains  were  able  to  relax  in  the  direction  of  the 

thermodynamic  state  of  higher  density. 

Both  series  of  fictive  temperatures  listed  in  Table  I  increased  with 

time  of  cure  t  and  achieved  values  greater  than  T  .  In  contrast,  the  Tr 

cure  6  c  f,a 

values  listed  in  Table  II  for  the  fully-cured  resin  decreased  with  time  of 

aging  t  .  at  130°C.  Figures  12  and  13  can  be  used  to  explain  these 
6  6  aging 

opposing  trends. 

•  Events  encountered  during  curing  at  a  temperature  T  <  T  (oo)  are 

$ 

depicted  in  Figure  12.  Curve  1  illustrates  a  hypothetical  cooling  curve 

(l°/min.)  obtained  for  the  reactants  without  reaction.  If  the  reaction  is 

chosen  to  begin  at  point  A,  T  >■  T  (the  T  of  the  inital  reactant 

c  g,r  g 

mixture).  After  substantial  reaction,  the  specific  volume  or  enthalpy  H 

of  the  system  decreases  from  point  A  to  B  and  the  material's  T^  becomes 

greater  than  Tc»  If  at  this  point  the  reaction  is  halted,  curve  2  would 

depict  the  resin's  heating  curve  without  accompanying  further  reaction.  Since 

during  part  of  the  curing  reaction  the  material  was  held  below  its  T  , 

8 

densif ication  (physical  aging)  occurs  simultaneously. 

During  the  first  heating  from  point  E  to  C  along  curve  2  the  densified 

glass  responds  sluggishly  and  endothermic  peaks  in  the  expansion  coefficient 

and  heat  capacity  are  observed.  The  maxima  of  these  peaks  are  at  the  point 

of  inflection  in  curve  2.  The  fictive  temperature  of  the  aged  glass  is  Tf 

r ,  a 

If  the  resin  is  cooled  from  point  C  to  D  before  appreciable  reaction  occurs, 

upon  immediate  reheating  at  a  rate  equal  to  the  previous  rate  of  cooling,  the 

T^  ^  is  an  approximate  measure  of  the  resin's  T^. 

As  the  curing  reaction  proceeds,  not  only  do  the  values  of  T^  and  T^  a 

increase,  but  the  difference  T,  -T,  also  increases  since  aging  increases. 

t  >  g  1 1  a 

This  is  manifested  by  the  increasing  area  of  the  attendant  endothermic  peak  in 
the  Cp-temperature  curve. 

The  decrease  in  the  T,  from  aging  without  attendant  chemical  reaction 

1 1  a 

is  illustrated  in  Figure  13.  Upon  cooling  from  temperatures  above  T  at 

o  ^ 

1  /minute  curve  3  is  generated.  If  the  cooling  was  halted  at  point  F 
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isothermal  volume  contraction  ensues.  When  the  specific  volume  v  or 

sp 

enthalpy  H  decreases  to  point  G,  the  T,  becomes  T,  ,  and  when  the  volume 

1 ,  a  i » i 

decreases  to  point  H,  a  lower  ^  *s  reached,  etc. 

The  characteristic  jump  in  the  heat  capacityAC^  occurred  at  the  . 

Although  not  apparent  in  Figure  4,  the  AC^  decreased  as  curing  progressed. 

The  AC  values  obtained  for  samples  of  various  cure  states  are  plotted  in 

Figure  14.  The  uncured  and  fully-cured  Epon  828/MDA  having  Ac^  values  of 

approximately  0.200  and  0.115  cal/gm.«°C,  respectively.  Bair  (16)  in  studying 

an  amine-cured  DGEBA  epoxy  found  the  AC  of  uncured  and  fully-cured  material 

to  decrease  from  0.123  to  0.089  cal/gm. •  C. 

The  measurements  of  Stevens  and  Richardson  (17)  show  that  the 

temperature  dependence  of  the  of  epoxy  glasses  is  greater  than  that  of  the 

liquid  of  rubber  above  T  .  As  curing  progresses,  with  the  concomitant 

8 

increase  in  T  ,  it  appears  that  the  entropy  of  the  glass  is  approaching  that 
8 

of  the  liquid/rubber.  The  decrease  in  AC^  appears  therefore  to  have  more 

dependence  upon  the  greater  rotational  and  vibrational  degrees  of  freedom  of 

the  glass,  and  not  the  loss  of  these  above  T  . 
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Graphical  Method  for  Determining  the  Fictive  Temperature  T* 

r ,  a 

The  fictive  temperature  T-  of  partially  and  fully-cured  material  can 

I ,  a 

be  approximated  graphically  from  endothermic  peak  areas.  Figure  15  shows  the 

first  DSC  specific  heat  scan  for  an  Epon  828/MDA  sample  cured  for  nine  hours 

at  130°C.  After  the  endothermic  peek  area  is  measured,  a  parallelogram  is 

constructed  to  the  left  of  the  peak.  Sides^)and(2) of  the  parallelogram  are 

parallel  to  th*’  curve  proceeding  the  peak.  Side(^,  the  right-side  boundary 

of  the  parallelogram,  is  located  on  the  onset  of  the  peak.  Side(4)will  be 

located  such  trat  the  area  of  the  parallelogram  equals  that  of  the  endothermic 

peak.  This  left-side  boundary  is  the  T^  ^  of  the  material. 

This  procedure  was  applied  to  the  specific  heat  scans  for  several 

partially  and  fully-cured  Epon  828/MDA  samples  aged  at  130°C  and  Dylene  8 

polystyrene  samples  aged  for  1,  2,  3,  6,  and  24  hours  at  85°C.  In  all 

instances  the  Tf  values  determined  by  this  graphical  method  were  within 
o  1  ,a 

1.5  C  of  those  determined  from  enthalpy  values. 
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CONCLUSION 


In  summary,  the  following  may  be  stated: 

1.  Differential  Scanning  Calorimetry  (DSC)  and  Infrared 

Spectrophotometry  (IR)  were  usee  to  monitor  the  degree  of  cure  within  Epon  828 

epoxy/MDA  samples  cured  at  130+^1°C.  Within  experimental  uncertainty,  the  data 

points  for  the  degree  of  cure  obtained  from  heats  of  reaction,  T,  and 

i » 6 

infrared  absorbance  as  a  function  of  time,  were  within  experimental  agreement 
with  one  another.  The  greatly  reduced  reaction  rate  during  the  final  20 
percent  of  cure  was  not  only  a  consequence  of  vitrification,  but  also  a  result 
of  the  depletion  in  the  number  of  reactive  epoxide  groups. 


2.  DSC  was  also  used  to  measure  the  extent  of  "physical  aging"  which 

took  place  during  the  cure  of  Epon  828/MDA  samples  at  130+1°C;  and  after, 

fully-cured  samples  were  aged  at  130°C,  37°C  below  their  ultimate  Tf 

*  >  8 


3.  Spontaneous  densification  occurred  whether  or  not  a  reaction  took 

place.  After  curing  for  one  hour  or  more  the  Tr  values  exceeded  the  T  a 

£>g  c 

a  dense  glass  was  formed.  The  glass  responded  sluggishly  to  temperature 
changes  near  the  fictive  temperature  and  upon  first  heating  an  endothermic 
peak  was  observed. 


'  A.  The  fictive  temperature  is  thought  to  be  the  hypothetical 

temperature  at  which  the  structure  would  be  in  equilibrium  and  was  determined 

from  DSC  heating  curves.  Two  sets  of  fictive  temperatures  were  obtained.  The 

values  of  T,  were  determined  from  enthalpy  data  obtained  during  first  DSC 

heating  scans.  They  reflect  the  amount  of  physical  aging  which  occurred 

during  the  curing  reaction,  after  the  material's  T^  was  greater  than  the 

The  correspondingly  higher  T,  values  were  obtained  directly  from  the  second 

* « 8 

heatings  after  aging  had  been  erased.  A  sample’s  T,  is  expected  to  be  close 

*  *  8 

to  its  T^  (normally  obtained  from  cooling)  if  the  heating  follows  immediately 
after  cooling  from  above  T  . 
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5.  The  endothermic  peak  temperatures  and  areas  for  the  epoxy  samples 

increased  with  curing  and  aging  time.  However,  enthalpy  data  revealed  that 

while  the  T,  of  partially-cured  samples  increased  with  curing,  the  T,  of 
I  ^  *  »3 

fully-cured  samples  decreased  with  aging  at  130  C. 

6.  The  reported  curing  and  aging  studies  revealed  that  although  longer 
periods  of  time  were  necessary,  fully-cured  samples  eventually  reached  a  level 
of  aging  equivalent  to  that  attained  by  samples  during  their  cure. 

7.  The  characteristic  jump  in  the  heat  capacity  AC  which  occurred  at 

the  T,  decreased  as  curing  progressed.  This  decrease  appears  to  be  dependent 
*  5  S 

upon  the  rotational  and  vibrational  degrees  of  freedom  of  the  glass. 

8.  A  graphical  method  of  determining  the  fictive  temperatures  Tf  for 

t ,  a 

partially  and  fully-cured  materials  from  measured  endothermic  peak  areas  was 
developed.  Values  determined  from  this  method  were  within  1.5°C  of  those 
obtained  from  enthalpy  data. 

9.  Physical  aging  is  assumed  here  to  be  the  consequence  of  the  material 

contracting  towards  its  equilibrium  density.  The  basic  idea  is  that  the  chain 

mobility  in  a  closely  packed  system  is  primarily  determined  by  the  degree  of 

packing  within  the  system.  Above  T  the  material  is  in  thermodynamic 

8 

equilibrium.  In  contrast  to  chemical  aging  or  degradation,  physical  aging  is 

a  thermo-reversible  process.  By  reheating  an  aged  material  above  its  T  the 

8 

original  thermodynamic  equilibrium  state  is  recovered  and  the  effects  of  aging 
erased.  The  effects  of  aging  can  then  be  repeatedly  reproduced. 
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Table  II 

Differential  Scanning  Calorimetry 
of  Physically  Aged  Epon  828/MDA  Samples 


Time  of  Aging 
log  sec. 


Peak 

Temp.,  °C 


Peak  Area 
Calories/Gram 


Fictive  Temp. ,  C 
Tf,a 


Approximate 

T-  a  (Note  1) 
t  ,  a _ 


4.76 

175.5 

0.20 

4.86 

176.8 

0.26 

5.54  • 

179.3 

0.55 

5.78  " 

182.0 

0.85 

6.24 

184.0 

1.15 

6.37 

.184.5 

1.25 

166.2 

167.5 

166.0 

164.3 

165.5 

163.5 

165.0 

163.0 

163.0 

162.4 

Note  1.  Graphical  "parallelogram"  method  used  to  determine  these  values. 
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Figure  4:  Fictive  Temperature  Determination  for  Epon  828/MSA  Samples  Cured  at  130+1°C. 

Second  DSC  Heating.  Heating  Hate  of  20°C/Min.  and  Cooling  Rate  at  80°C/Min. 

Sensitivity  of  0.5  mcal/sec.  full  scale.  Vertical  shift  arbitrarily  applied. 
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Figure  5.  Change  in  enthalpy  AH  versus  temperature  for  a  Epon  828/MDA 
samDle  cured  for  9. hours  at  130°C.  First  DSC  heating. 
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Figure  10.  Endothermic  peak  temperatures  or  Epon  828/MgA  samples 
after  curing  O  and  physical  aging  □  at  130  C. 
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Figure  11.  Endothermic  peek  arees  of  Epon  828/MDA  sampV s  after 
curing  Q#nd  physical  aging  □  at  130  °C. 


Figure  12.  Explanation  of  the  increase  in  and  for 
partially-cured  material  at  T  .  ’ 


Figure  13.  Explanation  of  the  decrease  In  Tf  for  fully-cured 
material  with  physical  aging  at  ’  T  *  T-. 
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Figure  15.  Graphical  method  of  determining  the  fictive  temperature  from  endothermic 
peak  area.  Epon  828/MDA  sample  cured  for  9  hours  at  ’  130°C.  First 

DSC  specific  heat  scan.  Heating  rate  of  20°C/min.  Sensitivity  of  0.5  mcal/sec 
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VOLUME  DEPENDENT  RATE  PROCESSES 
IN  AN  EPOXY  RESIN 

ABSTRACT 

Results  for  four  different  volume  dependent  rate  processes  of  an  epoxy  resin,  Epon  100 IF, 
fully  cured  with  a  stoichiometric  amount  of  4,4’  diamino  diphenyl  sulfone  are  presented:  1)  specific 
volume  measurements  at  constant  rates  of  cooling;  2)  time  dependent  volume  measurements  after 
rapid  temperature  changes;  3)  elongational  creep  compliance  measurements;  and  4)  shear  creep 
compliance  measurements.  Voluminal  and  shear  retardation  spectra  are  compared  as  are  the 
temperature  time  scale  shift  factors,  ar 
INTRODUCTION 

Kinetic  processes  in  amorphous  materials  are  linked  to  the  glass  temperature  because  Tg  is 
a  reflection  of  the  local  mobility  of  the  molecules  or  molecular  segments.  It  is  properly  determined 
in  a  cooling  experiment  because  it  is  close  to  the  first  departure  from  an  equilibrium  state  of  the 
liquid  being  cooled;  that  is,  at  higher  temperatures  as  soon  as  the  cooling  is  stopped  and  a  uniform 
temperature  is  achieved  the  random  fluctuating  disposition  of  the  molecules  is  fixed.  Every  polymer 
scientist  we  know,  including  ourselves,  is  aware  of  the  fact  that  Tg  is  a  function  of  the  rate  of 
cooling,  but  cannot  identify  a  set  of  cooling  experiments  that  establishes  this  fact.  This  behavior,  of 
course,  can  be  inferred  from  Kovac’s  isochronal,  0.02  hr.  and  100  hour,  glass  'rolume  -  temperature 
lines,  which  were  obtained  from  a  series  of  temperature-jump  experiments1.  The  polyvinylacetate 
studied  was  quenched  from  equilibrium  above  Tg  to  different  temperatures  between  -20°  and  35°C. 
It  is  not  possible  to  move  directly  from  point  to  point  along  the  isochronal  glass  line. 

It  was  desired  in  this  investigation  to  obtain  a  set  of  cooling  curves  for  a  fully  cured  epoxy 
resin  and  to  obtain  data  on  the  kinetics  of  isothermal  contraction  at  temperatures  where  the  specimen 
would  be  at  an  equilibrium  volume  before  the  rapid  change  in  temperature  was  made  to  the 
temperature  where  the  time-dependent  contraction  or  expansion  was  to  be  observed.  Since  the  rate 
of  change  in  volume  observed  near  Tg  is  predominantly  a  function  of  the  volume  itself  (specifically 
the  relative  or  fractional  free  volume),  the  measurements  of  such  changes  are  inherently  non-linear. 
This  has  been  made  abundantly  clear  by  the  classic  work  of  Kovacs1  and  others2,3.  The  non-linearity 
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was  reduced  in  this  study  by  keeping  the  temperature  change  at  2.5 °C.  The  isothermal  volume 
change  in  all  but  one  experiment  was  followed  until  the  equilibrium  density  was  reached. 

In  addition  to  dilatometric  measurements,  shear  and  elongational  creep  behavior  of  the  same 
glycidyl  ether  based  epoxy  resin  were  determined  over  a  range  of  temperatures.  Temperature  shift 
factors  were  obtained  for  the  four  different  kinds  of  measurements. 
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MATERIAL  AND  INSTRUMENTATION 

The  shell  epoxy  resin,  EPON  1001F  cured  with  a  stoichiometric  amount  of  4,4’  diamino 
diphenyl  sulfone  has  been  characterized  at  the  University  of  Akron  by  LeMay,  Swetlin  and  Kelly4. 

Tensile  creep  and  recovery  measurements  which  yielded  the  recoverable  tensile  creep 
compliance,  D(t),  above  and  below  Tg  were  carried  out  with  a  differential  elongational  creep 
apparatus5. 

1 

i  Torsional  creep  and  recovery  measurements  which  yielded  the  recoverable  torsional  creep 

I 

!  compliance,  Jp(t),  above  Tg  were  carried  out  with  a  frictionless  magnetic  bearing  torsional  creep 

j  apparatus  (Time-Temperature  Instruments,  Inc.,  P.O.  Box  40156,  Pittsburgh,  PA  15201). 

Dilatometric  measurements  were  made  using  a  modified  version  of  a  mercury  dilatometer. 

!  The  dilatometer  consisted  of  28  cm  of  0.05  cm  diameter  precision  bore  capillary  tubing  connected 

to  a  thin  walled,  stainless  steel  cylindrical  sample  chamber  (ID  =  1.5  cm,  height  =  4.0  cm,  wall 
thickness  =  0.031  cm).  The  sample  chamber  is  sealed  at  the  top  via  compression  of  a  teflon  O-Ring 
by  a  threaded  stainless  steel  cap.  Unlike  conventional  glass  or  quartz  diiatometers,  this  design  allows 
for  easy  installation  and  removal  of  a  sample.  The  glass  capillary  is  graduated  every  0.1  cm  allowing 
for  mercury  height  measurements  to  +  0.01  cm.  At  the  top  of  the  glass  capillary  there  is  a  ground 
glass  joint  for  attachment  to  the  mercury  distillation  apparatus.  Mercury  is  distilled  into  the 
dilatometer  at  a  pressure  of  less  than  5  microns  of  Hg. 

Phys;cal  aging  experiments  were  performed  on  the  EPON  1001F  near  its  glass  transition 
temperature  in  which  the  EPON  100 IF  sample  experienced  a  "fast"  temperature  jump  of  2.5°C  in  the 
temperature  range  124.0  -  136.5°C.  This  was  done  with  two  thermostatted,  well  stirred  silicon  oil 
baths  differing  in  temperature  by  2.5°C.  The  dilatometer  was  placed  in  one  bath  allowing  the  sample 
to  attain  its  thermal  equilibrium  volume.  The  dilatometer  was  then  transferred  to  the  other  silicon 
oil  bath  and  the  mercury  column  height  was  monitored  with  time. 

When  first  performing  these  aging  experiments,  it  was  observed  that  the  EPON  1001F  took 
an  unusually  long  period  of  time  to  thermally  equilibrate  to  the  bath  temperature.  This  was  due  to 
the  fact  inat  the  EPON  1001 F  sample  was  a  solid  cylinder  (length  =  3.8  cm,  Dia  =  1.30  cm)  with  a 
low  thermal  diffusivity.  EPON  100 IF  sample  was  removed  from  the  dilatometer  and  1/16"  diameter 
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holes  were  drilled  perpendicular  to  one  another  along  the  entire  length  of  the  sample.  This  allowed 
mercury  to  come  in  contact  with  the  sample  interior  increasing  the  heat  transfer  and  thus  decreasing 
the  time  for  the  sample  to  thermally  equilibrate.  With  this  "Swiss  cheese"  sample,  it  was 
experimentally  determined  that  the  EPON  1001F  sample  required  approximately  100s  to  thermally 
equilibrate  during  a  temperature  jump  of  -2.5*C  (T-Jump  down).  The  EPON  1001F  sample  required 
about  200s  to  equilibrate  during  a  temperature  jump  of  +  2.5'C  (T-Jump  Up).  This  significant 
increase  in  time  was  due  to  a  cooling  effect  from  ambient  air  that  occurred  in  transferring  the 
dilatometer  from  the  colder  bath.  During  T-Jump  Down  experiments  this  cooling  from  ambient  air 
aided  in  lowering  the  sample  temperature  thus  decreasing  the  time  for  the  sample  to  become 
isothermal. 

Specific  volume  of  the  EPON  100 IF  as  a  function  of  temperature  was  obtained  from  cooling 
experiments  over  a  temperature  range  of  T  -20  <  T  <  T  +  10°C.  The  sample  was  heated  to  T  + 

•*  9  9 

10°C  and  allowed  to  attain  thermal  equilibrium.  Once  at  equilibrium,  the  EPON  100 IF  was  cooled 
at  a  constant  rate  to  Tg  -20°C.  Time,  temperature,  and  mercury  column  height  were  monitored  to 
assure  that  the  cooling  rate  remained  constant. 

Temperature  control  was  achieved  by  modification  of  Bailey  constant  set-point  temperature 
controller.  A  motor,  time  relay  and  gearbox  were  added  to  decrease  the  set  point  at  a  constant  rate. 
The  gear  box  (Insco,  Main  Street,  Graton,  MA  01450)  possessed  turn  ratios  of  1/1  to  1/1000  to 
achieve  cooling  at  the  lower  rates.  A  timing  relay  was  added  to  turn  the  motor  on  for  a  set 
percentage  of  a  given  amount  of  time.  With  these  modifications,  constant  cooling  rates  covering 
almost  3  orders  of  magnitude  were  achieved. 

RESULTS 

Cooling  ai  Different  Rate? 

The  Epon  1001F  epoxy  resin,  fully  cured  with  DDS,  was  cooled  in  a  dilatometer  over  a  three- 
hundred-fold  range  of  rates;  from  0.90'  to  0.003'C  per  min.  Cooling  was  initiated  in  each  case  after 
equilibration  at  ~  140'C.  The  specific  volume  results  obtained  are  shown  in  Fig.  1.  The  equilibrium 
specific  volume  line  is 
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V(cm3/g)=  0.85446  +  4.28  X  10'4  (T  -  130°C)  (1) 

The  liquid  -  like  and  glassy  cubical  thermal  expansion  coefficients  at  Tg  are  at  (Tg)  =  5.01  X  10’4 
and  ag  (Tg)  =  1.93  X  10'4/°C.  Tg  is  identified  with  the  intersection  of  the  equilibrium  and  glassy 
lines.  The  specific  volumes  at  the  different  rates  of  cooling  Q  =  3.0  X  10'3,  5.0  X  10'2,  2.5  X  10'1, 
and  9.0  X  10'1oC/min  at  Tg(Q)  of  126.1%  129.9",  130.5",  and  131.7"C  are  0.85277,  0.85402,  0.85469 
and  0.85523  cm3/g,  respectively.  The  results  of  V(Tg,Q)  are  probably  relatively  accurate  to  several 
parts  per  1 0s  and  absolutely  to  about  one  part  in  104. 

The  temperature  range  over  which  measured  specific  heat  values  depart  from  equilibrium  and 
glassy  lines  in  the  absence  of  physical  aging  (densification)  effects  has  been  called  the  "glass 
transformation  range".  This  temperature  region  corresponds  to  that  where  measured  specific  volume 
values  lie  above  the  equilibrium  and  glassy  lines.  This  range  of  departure  is  a  reflection  of  the 
breadth  of  the  viscoelastic  volumetric  retardation  function6,7.  As  such  it  should  be  a  function  of  the 
rate  of  cooling.  Careful  examination  of  Fig.  1  shows  this  to  be  so.  If  the  departures  of  the  data 
points  from  the  lines  (AV  =  Vm#as_  -Vl{nes)  are  plotted  as  a  function  of  temperature,  it  is  eminently 
clear  that  the  "transformation  range"  decreases  appreciably  with  a  decreasing  rate  of  cooling;  see  Fig. 
2.  The  half-width  decreases  by  about  a  factor  of  two  from  the  highest  to  the  lowest  rate.  As  the  rate 
decreases,  mechanisms  of  contraction  with  shorter  retardation  times  associated  with  the  changing 
liquid  structure  have  time  to  reach  their  equilibrium  state.  At  the  lower  temperatures  all  of  the 
cooling  curves  reach  the  glassy  lines  at  about  the  same  time  because  the  temperature  shift  factors  aT 
of  the  mechanisms  are  growing  at  a  greatly  accelerating  rate  with  decreasing  temperature  because  of 
the  rapidly  increasing  molecular  crowding. 

Temperature  Jump  Experiments 

When  the  temperature  of  an  amorphous  material  is  increased  or  decreased  rapidly  near  or 
below  Tg,  after  thermal  equilibrium  is  established  the  volume  of  the  material  normally  continues  to 
change  in  the  direction  of  the  equilibrium  value  for  the  new  temperature.  Near  and  below  Tg  some 
of  the  molecular  motions  involved  in  attaining  the  equilibrium  structure  are  very  sluggish  because 
of  their  complexity  in  the  crowded  state.  These  are  believed  to  be  the  same  molecular  motions  that 
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are  involved  in  the  approach  to  an  equilibrium  liquid  state  following  a  change  in  the  pressure  in  the 
same  temperature  range.  Rehage  and  Goldbach  have  shown  the  similar  response  of  a  polystyrene  to 
temperature  and  pressure  jumps3,8.  The  response  to  temperature  changes  should  be  recognized  to 
be  a  volumetric  viscoelastic  response. 

Kovacs1  has  investigated  time  dependent  volume  changes  following  several  different  kinds 
of  thermal  histories: 

1)  temperature  jumps  from  equilibrium  above  Tg  at  a  specific  temperature  to  various 
temperatures  below  Tg; 

2)  temperature  jumps  from  equilibrium  at  a  specific  temperature  below  T  to  several 
temperatures  above  Tg; 

3)  jumps  from  equilibrium  at  various  temperatures  above  Tg  down  to  the  same 
temperature  below; 

4)  jumps  from  equilibrium  at  several  temperatures  below  T  up  to  a  specific  temperature 

above  Tg;  ' 

5)  multiple  temperature  jumps  with  memory  reversals  of  the  volume  changes. 

The  temperature  jump  experiments  that  are  being  reported  here  are  variations  on  Kovac’s 
histories  1)  and  2).  In  our  experiments  the  Epon  1001F  was  either  cooled  or  heated  by  2.5°C  after 
the  fully  cured  epoxy  reached  its  equilibrium  density  and  the  change  in  volume  was  observed  until 
equilibrium  was  achieved  before  an  additional  step  was  taken.  The  sole  exception  was  the  last  cooling 
step,  where  the  time  to  reach  an  equilibrium  density  was  excessive. 

An  example  of  a  pair  of  curves  expanding  and  contracting  toward  equilibrium,  following 
2.5°C  temperature  steps  is  shown  in  Fig.  3.  The  volume  changes  following  temperature  shifts  from 
129*  and  134°  to  131.5*C  are  shown.  The  ordinate  is  actually  the  fractional  departure  from  the 
equilibrium  specific  volume  curve.  Both  curves  take  the  better  part  of  a  day  to  reach  equilibrium. 
However,  two  differences  can  be  seen  between  the  contraction  and  expansion  curves.  The  asymmetry 
can  be  seen  which  reflects  the  intrinsic  non-linearity  of  these  experiments.  The  ascending  curve  is 
delayed  in  response  because  of  the  great  molecular  crowding,  but  reaches  greater  absolute  slopes 
because  of  the  autocatalytic  nature  of  the  expansion.  The  contraction  curve  descends  at  an  earlier 
time  toward  equilibrium,  but  decelerates  relatively  more  because  a  state  of  greater  crowding  is  being 
approached.  The  time  of  contraction  and  expansion  are  corrected  for  the  approximate  time  to 
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achieve  an  equilibrium  temperature  after  the  dilatometer  is  transferred  from  one  thermostatted  bath 
to  the  other;  i.e.  tf  is  the  approximate  time  for  the  temperature  to  become  uniform  throughout  the 
dilatometer  bulb  after  the  transfer.  Kovacs  was  able  to  assume  a  t,-  =  100  sec.  for  both  cooling  and 
heating  experiments  because  his  transfers  of  polyvinylacetate  dilatometers  were  made  near  room 
temperature.  We  were  required  to  use  a  larger  thermal  equilibration  time  for  the  transfers  to  the 
higher  temperature  bath,  since  at  elevated  temperatures  ~  130°C,  cooling  occurred  during  the  several 
seconds  required  for  the  transfer.  Thermal  equilibration  was  thus  delayed,  whereas  during  the  step 
down  transfer  the  equilibration  time  was  diminished  by  the  extra  cooling  incurred  during  the 
transfer.  The  tt  values  were  determined  at  temperatures  somewhat  above  Tg  where  the  epoxy  resin’s 
retardation  times  were  a  few  seconds  or  less. 

The  absolute  fractional  departures  from  equilibrium  obtained  at  four  temperatures  after 
heating  from  equilibrium  and  five  temperatures  after  cooling  from  equilibrium  are  presented  in  Fig. 

4.  Two  and  a  half  degree  temperature  steps  were  used  in  every  case. 

With  134.0°C  chosen  as  the  reference  temperatures  both  sets  of  data  were  reduced  to  common 
curves  which  range  over  an  extended  logarithmic  time-scale.  The  reduced  curves  are  shown  in  Fig. 

5.  The  reduced  curve  obtained  following  cooling  steps  is  more  extended  simply  because  of  the  long 
contraction  times  involved  at  124*C.  At  126.5*C  it  took  about  3  to  4  days  to  reach  equilibrium. 
Therefore,  a  corresponding  heating  curve  starting  at  124°C  could  not  be  obtained.  Only  the  reduced 
contractions  curve  will  be  further  analyzed  here.  The  very  same  sample  and  dilatometer  were  used 
in  both  the  constant  rate  of  cooling  and  temperature  jump  measurements. 

Determination  fif  jh£  Bulk  Retardation  Spectrum 

The  reduced  curve  of  the  fractional  volume  departure  from  equilibrium  which  was  obtained 
following  cooling  steps  [V(t)  -  V(oo)]/V(oo)  was  normalized  with  the  factor  (Voo/[V(0+)  -  V(oo)]  yielding 
the  fraction  of  the  volume  change  at  time  t:  [ V(t)  -  V(oo)]/[V(0*)  -  V(oo)].  V(t)  is  the  specific  volume 
at  t;  V(oo)  is  the  equilibrium  specific  volume  and  V(0+)  is  the  volume  obtained  following  an 
instantaneous  temperature  jump.  It  is  related  to  the  equilibrium  volume  before  the  jump  V(O)  by 

V(0*)  =  V(0)(l-a9  AT),  (1) 

where  ag  is  the  thermal  expansion  coefficient,  3(ln  V)/3T,  of  the  glass. 
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Volume  changes  can  be  caused  by  changes  in  temperature  or  pressure.  The  changes  in  volume 
are  related  to  the  compressibility  by7 


|[V(0+)  -  V(oo)]/[V(0)  -  V(oo)U[\  Bit.). ;  B  1  V(t)  -..V£oo) 

1  JL1_B(oo)-U  -  V(0)-  V(oo) 


(2) 


where  B(t)  is  the  time  dependent  compressibility,  -  3[tn  V(t)]/3P.  B(oo)  and  Bg  are  the  equilibrium 
and  glassy  compressibilities,  respectively.  For  pressure  changes  [V(0+)  -  V(oo)]/[V(0)  -  V(oo)]  =  [B(oo) 
-  Bg]/B(oo)  and  for  temperature  steps  [a(oo)  -  ag]/a(oo). 


For  our  reduced  contraction  curve. 

,  V(tT-  V(oo).  _  B(t)-B 

'  V(0+)  -  V(ooj  "  BI^pBg 

or  if  we  define  a  normalized  volume  retardation  spectrum  Ly(r)*  where 
f  L  (r)  d  In  r  =  I 


(3) 


(4) 


then 


YLQ*)  -  YCtl 
V(0*)  -  V(oo) 


=  JVv  (f)(l-e’t/T)  d  In  t. 

-oo 


(5) 


Using  iterative  computer  calculations9  Lv(r)  has  been  obtained  from  the  normalized  reduced  data. 
The  resuPant  bulk  spectrum  is  shown  in  Fig.  6  where  it  is  compared  with  the  shear  spectrum  which 
was  deduced  from  the  reduced  shear  creep  compliance  J(t)  curve10.  The  abscissa  for  log  Lv  (r)  was 
chosen  so  that  the  short  time  portions  of  the  spectra  superpose.  The  fact  that  they  can  be 
superimposed  for  at  least  two  decades  of  time  scale  suggest  that  the  same  molecular  motions  are 
involved  in  both  the  shear  deformation  and  the  densification.  The  longer  time  mechanisms  which 
involve  larger  range  coordinated  polymer  chain  backbone  motions  are  reflected  in  the  long-time  peak 
in  the  shear  spectrum.  It  is  no  surprise  that  the  long  range  molecular  motions  have  no  measurable 
effect  on  the  local  ordering  that  determines  a  material’s  density. 


L v(r)  d  In  r  is  the  normalized  distribution  f unction  of  volume  retardation  times  lying  between 
In  r  and  In  r  +  d  In  r  which  contribute  to  the  fraction  of  the  equilibrium  bulk  (voluminal) 
compliance,  B(oo)  -  Bg. 
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Creep  Measurements 

Elongational  creep  compliance,  D(t),  measurements  were  carried  out  with  a  differential 
instrument5  in  which  the  lengths  of  twin  samples  were  monitored.  One  sample  was  subjected  to  the 
tensile  creep  loading  and  fluctuations  of  length  due  to  thermal  variations  cancelled  out  in  the 
differential  signal.  Reduced  Log  Dr(t)  curves  obtained  at  10  temperatures  between  117.7*  and 
150.2*C  are  presented  in  Fig.  7.  The  small  rubberlike  temperature  correction  was  applied  to  the  time 
dependent  contribution  above  the  glassy  compliance, 

Dp(t)  -  [DU)  -  DgllTp/VJ  +  D9.  (6) 

T  and  T0  are  the  absolute  temperatures  of  measurement  and  reference,  130°C,  respectively,  p  and  p0 
are  the  corresponding  densities.  The  corrected  curves  were  reduced  by  horizontal  time-scale  shifts 
onto  the  curve  measured  at  130*C.  At  this  chosen  reference  temperature  the  material  had  its 
equilibrium  density.  To  reach  equilibrium  a  day  was  required  as  indicated  by  the  above  contraction 
curves.  Excellent  superposition  of  the  curve  was  achieved  as  can  be  seen  in  Fig.  8.  The  Dp(t) 

reduced  curve  is  compared  with  the  reduced  shear  compliance  Jp(t)  curve  shown  as  a  long-dashed 

line.  The  torsional  creep  measurements  from  which  the  Jp(t)  curve  was  obtained  are  new 
measurements  which  are  in  close  agreement  with  previous  measurements  on  fully  cured  EPON 
1001F/DDS.10  The  same  mix  was  used  in  preparing  the  specimens  for  both  the  torsional  and 
elongational  measurements  in  hope  that  some  information  on  the  bulk  compliance  B(t)  could  be 
obtained  from  the  relation11 

B(t)  *  9  D(t)  -  3J(t)  (7) 

at  least  at  relatively  short  times  where  D(t)  and  J(t)  are  not  much  larger  than  B(t).  Comparison  at  long 
reduced  times  >  104  sec.  where  the  influence  of  B(t)  is  not  discernible  indicated  that  the  Jp(t)  curve 
was  at  longer  times  by  0.8  of  a  decade  of  time.  This  obviously  reflects  a  higher  glass  temperature  of 
the  specimen  on  which  the  torsional  measurements  were  made.  Beyond  reduced  times  of  103sec  Jp(t) 
and  Dp(t)  are  superposable  within  experimental  scatter  with 

Jp(t)  r  3Dp(t)  (8) 

indicating  a  Poisson’s  ratio  of  1/2.  Since  direct  measurements  by  I.  -C.  Chay  indicate  that  the 
equilibrium  bulk  compliance  B^  3  X  10'11  cmz/dyne  for  this  material,  no  detectible  dilation  is 
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expected.  The  long-dashed  line  for  J(t)  has  been  shifted  to  shorter  times  by  the  required  Alog  t  of 
0.80.  The  short-dashed  line  represents  log  Jr(t)  shifted  vertically  by  0.477,  the  logarithm  of  3. 
Values  for  B(t)  were  estimated  in  the  region  where  the  solid  and  short  dashed-lines  diverge.  They 
were  found  to  be  substantially  greater  than  the  directly  measured  equilibrium  value  by  about  50%. 
We  believe  this  is  a  result  of  distortions  that  occur  in  the  time-temperature  reduction  process. 
Extrapolating  both  curves  to  zero  time  using  a  fit  to  Andrade  Creep10  at  short  times,  a  glassy 
compliance  of  Dg  =  4.90  X  10'11  cm2/dyne  and  Jfl  =  1.36  X  10'10  cm2/dyne  are  obtained.  These  yield 
a  glassy  bulk  compliance  Bg  of  3.3  X  10'11  cm2/dyne  and  a  glassy  Poisson’s  ratio  of  0.35  from 

"g“(V2D9>-l  (9) 

These  values  are  both  clearly  too  large  for  reasons  unknown  to  us.  The  differences  involved  do 
magnify  experimental  uncertainties  by  an  order  of  magnitude. 

Discussion 

To  check  the  hypothesis  that  local  coordinated  molecular  motions  are  the  primary  contributors 
to  the  change  in  structural  arrangements  in  amorphous  materials  which  is  involved  in  time-dependent 
dilatation  or  densification,  the  kinetic  volume  changes  in  the  glass  transformation  range  have  been 
measured  and  the  results  have  been  compared  with  the  creep  compliance  curves  which  reflects 
coordinated  molecular  motions  over  a  very  wide  range  of  the  time-scale.  If  the  same  molecular 
motions  are  indeed  involved  in  density  changes  and  various  mechanical  deformations  the  temperature 
dependencies  must  be  the  same. 

The  changes  with  temperature  in  the  neighborhood  of  Tg  are  known  to  reflect  the  intrinsic 
dominating  variable,  the  free  volume.  Changes  in  rate  under  isothermal  and  isobaric  conditions  are 
testing  to  this  fact12'14.  The  comparison  of  temperature  shift  factors  shown  in  Fig.  9  does  indicate 
that  the  changes  in  volume  have  the  same  temperature  dependence  as  the  creep  compliances. 
Markovitz  has  shown11  that  reasoning  from  Eq.  7  that  for  the  three  compliance  functions  J(t),  D(t), 
and  B(t)  of  a  material  to  be  temperature  superposable  they  must  all  have  the  same  time-scale  and 
magnitude  shift  factors.  The  data  presented  above  may  be  the  first  such  test  and  are  certainly  in 
accord  with  his  conclusions. 
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In  addition  the  normalized  voluminal  retardation  spectrum  Ly(r)  calculated  from  the  reduced 
isothermal  contraction  curve  appears  to  correspond  to  that  obtained  from  the  shear  compliance  curve 
Ls(r)  at  short  times.  The  fact  that  Ly(r)  does  not  extend  into  the  realm  of  long-times  supports  the 
expectation  that  long-range  coordinated  molecular  motions  do  not  play  a  significant  role  in 
determining  the  liquid  structure  which  determines  the  density. 
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FIGURE  LEGENDS 

Fig.  1)  The  specific  volume  \{cm3/g)  of  EPON  1001F  fully  cured  with  a  stoichiometric  amount  of 
4,4/  -  diamino  diphenyl  sulfone,  DDS  shown  as  a  function  of  temperature  at  four  different 
rates  of  cooling  0.90,  0.25,  0.050,  and  0.003  °C/min.  Glass  temperatures  identified  by  the 
intersection  point  of  the  equilibrium  and  glass  lines  are  listed. 

Fig.2)  Deviation  of  measured  specific  volume  points  from  the  equilibrium  and  glass  lines  of  Fig.  1 
plotted  as  a  function  of  temperature,  showing  the  extent  of  the  '’transformation  range"  and 
its  change  with  the  rate  of  cooling. 

Fig.  3)  Fractional  volume  deviation  from  equilibrium  as  a  function  of  time  at  131.0°C  after  cooling 
from  134°C  and  after  heating  from  129°C  plotted  versus  the  logarithm  of  the  corrected  aging 
time  t-tj  where  tj  is  the  estimated  time  for  the  specimen  to  reach  a  uniform  temperature. 

Fig.  4)  Absolute  fractional  volume  deviation  from  equilibrium  as  a  function  of  the  logarithm  of  the 
corrected  aging  time  for  five  temperatures  after  cooling  2.5°C  from  equilibrium  and  for  four 
temperatures  after  heating. 

Fig.  5)  Volume  contraction  and  expansion  curves  from  Fig.  4  reduced  by  time-scale  shifts  to  the 
reference  temperature  T0  =  134°C.  The  time  t  is  the  corrected  value. 

Fig.  6)  Comparison  of  retardation  spectra  for  voluminal  Ly  and  shear  deformation  Ls.  In  this  double 
logarithmic  plot  of  the  distribution  functions  of  retardation  times  r  the  ordinate  scales  have 
been  adjusted  to  superpose  the  short  time  results. 

Fig.  7)  Double  logarithmic  presentation  of  the  reduced  elongational  creep  compliance  Dp(t)  as  a 
function  of  time  for  EPON  1001F  at  10  temperatures  as  indicated.  The  temperature  of 
reduction  is  130'C. 
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Fig.  8)  Elongational  creep  compliance  Dp(t)  curves  from  Fig.  7  shifted  to  superpose  at  a  reference 
temperature  T0  of  130°C.  The  reduced  shear  creep  compliance  Jp(t)  curve  also  shown  for 
comparison  (long-dashed  line).  The  latter  has  been  shifted  toward  shorter  times  by  0.8  of 
a  decade  to  match  long-time  rates.  The  short-dashed  line  is  the  Jp(t)  curve  diminished  by  a 
factor  of  three. 

Fig.  9)  Logarithm  of  time-scale  shift  factors  of  EPON  100IF  for  four  different  kinds  of  rate  process 
experiments;  1)  volume  -  constant  rate  of  cooling;  2)  volume  -  temperature  jump;  3)  shear 
creep  compliance  and;  4)  elongational  creep  compliance,  presented  as  a  function  of 
temperature. 
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Summary 

The  behavior  of  polymer-based  material  in  the  vicinity  of  a  crack  tip  and 
the  subsequent  growth  of  a  crack  through  the  material  are  topics  of  current 
research.  The  model  described  in  this  paper  is  based  on  a  finite  element 
analysis  of  the  region  surrounding  a  crack  tip.  The  parameters  used  in  the 
model  allow  for  a  variety  of  different  materials  including  those  displaying 
non-Hookian  behavior.  The  model  includes  a  simple  mechanism  for  examining 
fatigue  leading  to  crack  propagation.  Further  investigation  of  fatigue  crack 
growth  can  be  done  using  the  results  of  the  present  model. 
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I.  Introduction 

The  fundamental  assumptions  of  classical  theoretical  analysis  of 
rubber-like  elasticity  are  basically  as  follows: 

1.  The  distribution  of  the  chain  dimensions  is  describable 
by  the  Gaussian  statistics  of  freely-jointed  chainsJ 

2.  Intermolecular  forces  are  insignificant. 

3.  Internal  energy  changes  during  deformation  are  non-existent. 

4.  Network  chains  deform  in  proportion  to  the  strain  on  the 
sample  (affine  deformation). 

These  assumptions  lead  to  an  expression  for  the  force  per  unit  unstrained 
cross-sectional  area2 

f  =  G(X  -  A'2)  (1.1) 

where 

G  =  NKT  (1.2) 

X  =  L/Lq  =  1  +  e  (1.3) 

In  an  attempt  to  treat  a  wider  variety  of  polymers  than  the  case  of  ideal 
rubber,  the  first  assumption  can  be  avoided  by  using  a  result  of  Treloar's 
.  for  non-Gaussian  networks  given  by2'*^ 


f  =  XG(X  -  X-2) 


(1.4) 
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The  first  five  terms  of  X  are 

X  =  1  +  An-’^A2  +  4A"1) 

+  Bn“2(5A4  +  8A  +  8A'2) 

(1.5) 

+  Cn‘3(35A6  +  60A  +  72  +  64A-3) 

+  Dn"4(630A8  +  1120A  +  1440X  +  1536A"1  +  1280A'4) 

where  A  =  35/25; 

B  =  297/6125; 

C  =  12312/2205000; 

D  =  12611 7/ [(693) (673750)] ;  and 
n  is  the  number  of  random  links. 

It  is  easy  to  see  that  expression  (1.4)  reduces  to  the  Gaussian  case, 

(1.1 ),  for  large  n. 

The  proposed  finite  element  approach  enables  the  use  of  Equation  (1.4) 
indirectly.  That  is,  using  the  derivative  with  respect  to  extension  as  an 
expression  for  the  modulus  within  each  element.  This  aids  in  modeling 
overall  non-linearity.  To  demonstrate  simple  fatigue  behavior,  the  number 
of  effective  crosslinks,  N  in  Equation  (1.2),  is  expressed  as  a  function 
of  the  amount  of  strain  and  the  length  of  time  that  the  material  is  subjected 
to  that  strain.  This  is  given  by 

N(6t,e)  *  N(t-6t)exp(-oc6t)  (1.6) 

where  a  is  an  adjustable  parameter  which  has  dimensions  of  frequency. 

The  details  of  the  finite  element  method3,4  are  outlined  on  the 


following  pages. 
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II.  Finite  Element  Method 

Shown  in  Figure  1  is  a  discretized  slab  representing  the  material  in 
the  vicinity  of  a  crack  tip.  All  of  the  elements  are  rectangular  for 
ease  in  calculation  and  visualization. 

Consider  the  i^h  element  shown  below. 


At  each  corner  (node)  there  are  two  degrees  of  freedom,  i.e.,  displacements 
in  the  x  and  y  directions.  These  displacements  are  described  by 

u(x,y)  =  Ao  +  A-iX  +  A?y  +  A3xy 

(II. 1) 

v(x,y)  =  A4  +  A5X  +  Agy  +  A?xy 

where  u(x,y)  is  displacement  in  the  x  direction  and  v(x,y)  is  displacement 
in  the  y  direction.  The  form  of  Equation  (II. 1)  insures  continuity  across 
element  boundaries.  Note  that  for  constant  x  the  displacements  vary 
linearly  with  y,  and  similarly  for  constant  y.  Equation  (II. 1)  can  be 
written  in  matrix  form  as 

{6}  ="{£=[M]{A}  (II. 2) 


[M]  = 

1  X 

y 

xy 

0 

0  0  0  ' 

0  0 

0 

0 

1 

x  y  xy 

{A}1  = 

<Ao 

Ai 

CVJ 

A3 

A4  A5  Ag  A7) 

where 
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This  formulation  can  be  extended  to  represent  all  of  the  nodes  in  the  ith 


element  as  shown  below. 
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{6}  *  [c]{A} 

The  vector  {A}  can  then  be  given  by 

{A}  =  [c]-1^} 

The  strain  within  a  single  element  is  defined  in  the  usual  fashion 


(II. 3) 


(II. 4) 
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where 

*  A]  +  A$y  =  0  +  A]  +  0  +  A^ 

=  Ag  +  A7X  =  Q  +  0  +  Ag  +  A7X 

*  Ag  +  A3X  =  0  +  0  +  A2  +  A3X 

=  A5  +  A7y  =  0  +  A5  +  0  +  A7y 

Again,  this  can  be  expressed  in  matrix  form  by 


or 


1  0  y 
0  0  0 
0  1  x 


0  0  0  .0 

0  0  1  x 

0  1  0  y 


(II. 5) 


{£}  «  [N]{A> 

Finally,  substituting  Equation  (II. 4)  for  {A}  into  the  above  equation 
gives 


{e}  = 


(II. 6) 
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The  "stiffness"  matrix  [K]  can  be  obtained  by  requiring  that  the 
potential  energy  be  minimized  with  respect  to  the  displacements.  The 
potential  energy  for  this  case  is  given  by 

7Tp  =  ;v  Icodv  (II. 7) 

where 


and 


{0}  = 


yy 


=  Ed] {e} 


xy  ' 


(II. 8) 


[d] 


1 

(1  -  v2) 


1 

V 

0 


v  0 

1  0 


(II. 9) 


Substituting  expressions  (II. 6)  and  (II. 8)  for  e  and  o  into  (II. 7)  gives  for 
the  potential  energy 


up  =j/{6}t[c-1]t[N]t[d][N][c-'1]{5}dv  (II.  10) 

By  equating  the  partial  derivatives  of  7rp  with  respect  to  the  displacements 
of  zero,  i.e.. 
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the  following  expression  for  the  stiffness  matrix  is  obtained 

[K]{<5}  =  0 


where 

[K]  =  /vCc“1]t[N]t[d][N][c"1]dv  (11.11) 

This  process  is  continued  for  each  element  in  the  mesh  shown  in  Figure  1. 
Finally,  the  overall  stiffness  matrix  is  formed  from  the  individual  element 
stiffness  matrices. 

Before  solving  the  set  of  simultaneous  equations  for  the  displacements, 
boundary  conditions  must  be  implemented.  The  boundary  conditions  used  in 
this  model  are  as  follows: 

1.  The  displacement  in  the  vertical  direction  at  the  top  of  the 
slab  is  equal  to  the  negative  vertical  displacement  at  the 
bottom  of  the  slab  and  both  displacements  are  constant  in  the 
horizontal  direction. 

2.  The  horizontal  displacements  along  the  top,  bottom  and  sides  is 
zero,  except,  of  course,  in  the  area  representing  the  crack  tip. 

3.  The  vertical  displacement  along  the  axis  of  symmetry  is  zero. 

The  first  of  these  conditions  is  simply  that  the  slab  is  being  stretched 
by  equal  amounts  in  opposite  directions  at  the  top  and  bottom.  Since  no 
loads  are  specified,  the  displacements  are  given  as  boundary  conditions. 

The  second  condition  amounts  to  stating  that  the  slab  under  consider¬ 
ation  is  far  removed  from  the  edges  of  the  material. 

The  third  condition  is  based  on  the  assumption  that  the  geometric  axis 
of  symmetry  will  remain  stationary  in  the  vertical  direction  due  to  the 


symmetry  of  the  applied  strain.  This  assumption  allows  the  solution  of  the 
upper  half  of  the  slab  to  be  duplicated  for  the  lower  half,  thus  reducing 
the  number  of  equations  by  half. 
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III.  Time  Evolution 

To  model  fatigue  behavior,  the  displacement  at  the  top  and  bottom 
boundaries  is  cyclic.  In  this  particular  case  the  displacement  is  given 
by 

d  =  da  -  d5  cos(2rfn5t)  (III.l) 

where  f  is  frequency; 

n  =  0,1 ,2,3.. .ny,  nj  is  the  total  number  of  time  steps  per 
cycl  e ; 

6t  is  the  length  of  each  time  step; 

da,  dfc  are  adjustable  parameters  for  the  amplitude  of  the 
displacement,  da  >  djj. 

The  finite  element  method,  as  previously  described,  is  applied  at 
each  time  step  6t.  The  function  for  the  number  of  effective  crosslinks. 
Equation  (1.6),  is  updated  and  the  modulus  within  each  element  is 
calculated.  Then,  as  the  sample  evolves,  its  characteristics  can  be 
examined  at  any  point  in  time.  Measurements  could  be  made  of  a  variety  of 
quantities  such  as  stress  in  the  vicinity  of  the  crack  tip,  stored  elastic 
energy,  or  the  change  in  stored  energy  with  respect  to  the  formation  of 
new  surfaces. 

The  criteria  for  the  formation  of  new  surfaces  is  presently  based  on 
the  amount  of  strain  at  the  crack  tip.  Propagation  of  the  crack  is 
constrained  to  the  axis  of  symmetry. 

Graphics,  which  are  displayed  at  every  time  step,  allow  visualization 
of  the  changing  sample.  Figures  2-5  show  a  sample  at  stages  just  before 
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and  after  tearing.  Figures  6-9  show  the  strain  contours  during  the  same 
time  period. 
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Conclusions  and  Suggestions  for  Further  Work 

The  model  described  in  this  paper  can  be  used  to  gain  valuable 
information  about  the  behavior  of  material  in  the  vicinity  of  a  crack 
tip,  particularly  material  in  transition,  i.e.,  changing  modulus.  By 
varying  the  parameters,  totally  different  characteristics  can  be  displayed 
both  visually  and  quantitatively. 

It  should  be  noted,  however,  that  the  present  model  is  a  starting 
point  or  basis  for  further  investigation.  It  would  be  interesting,  for 
example,  to  introduce  different  fatigue  mechanisms  such  as  local  heating, 
or  tc  use  more  complicated  criteria  for  crack  propagation.  Ultimately, 
the  description  of  two-dimensional  crack  growth  should  be  attainable 
using  this  approach. 
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FIGURE  CAPTIONS 

Fig,  1  Discretized  slab  in  the  initial  state,  i.e.,  no  strain. 

There  are  264  total  elements.  The  area  which  has  been  removed 
represents  the  initial  crack  tip. 

Figs.  2-5  Vertical  strain  is  applied  in  equal  increments  for  each 

successive  figure.  Figures  4  and  5  show  the  formation  of  new 
surfaces. 

Figs.  6-9  Hand-drawn  regions  of  strain  intensity  based  on  the  strain  at 
element  boundaries  are  added,  to  Figures  2-5.  The  regions  are 
not  established  by  absolute  strain  but  rather  the  ratio  of  the 
strain  at  a  point  to  the  strain  at  the  crack  tip  on  the  axis 
of  symmetry,  where  the  maximum  strain  occurs. 

Figures  1-9  follow. 
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APPENDIX 

Source  Code  of  FORTRAN  Program 


c  main 

dimension  x  <200. 4 > , y  <200, 4) , s <200, B, 8),c<8,8),ct(8,8>,ck<3,8,8 
di  mension  p  <350  )  ,  del  <350)  ,  nl  <200, 4 )  ,  n2  <4 )  ,  hh  <  10)  ,  w  (20) 
real  nu.n  <8.  8) , k  <350,350) ,nec  <200) 
integer  g,bc<ll) 

3-0 

-flg-0 

t::-0 

read  <* , *> nts, dt , 4 , nu, h, nel , nei , nd-f , dx , da, db, dcr it , g 
call  inputl<hh,w> 
call  gennod (x,y,nl,hh, w) 
do  100  i* 1,200 

nec  < i > -1 . 0  e  23 

100  continue 
nen=161 

do  101  i =1.350 

101  del <i ) =6 
cal  1  mai n (  ) 

call  graphics  <:s ,  y,  nl ,  nel ,  del ,  tx  > 

1  call  cal cdi sp (nts, dt, f , disp, g, da. db ) 
do  5  i-1,350 

do  5  j-1,350 

5  k  < i , j  > -0. 

do.  Jt0..J.t*l,  nel _ _ _ I _ _ _ _ _  .  .. _ - 

call  calcnec  (nec, del ,  dt , nl , y ,  it)  . 

'  call  -forme  (del ,  e,  nu,  n  1,  y, nec*  it ) 

...  call  -formcinv  (c, ,  y,  ct,  it,  ck) 

_  call  -formn  (n, x  ,  y,  nu,  e,  h,  i  t> 

call  mmul t <ct , n, c, s, it ) 

call  sti-fmat  (k,s,nl,n2,it,  nen, nel ,  nd-f  , nei  > 

10  continue 

call  tstbc (be, del , y, nl , dcri t ) 
call  boundcon <k , di sp , p, bw. be) 

.  call  matsol v  <k, nen, p,  ndf ,  del ) 
call  mainO 

call  graphics ,y,nl , nel ,  del ,  tv.  > 
g=Q+i 

-  i-f  <g.  It. nts)  go  to  1  . 

g=0 

go  to  1 

stop 

end 

c  ' 

;  c--~  •  ■  v.’ . '  .... 

e-  - 

subroutine  ealedisp <nts,dt,f ,disp,g,da,db> 

integer  g 

pi-3.1413927 

di  sp-da-db*cos  (2*pi  *-f*dt*g> 
return 
end 
c 
e 

subroutine  tstbc  <bc,del,y,nl,dcrit> 
dimension  del (350) , y (200,4) ,nl (200,4) 
integer  bc<ll) 
do  40  i-1, 10 
j-31«-(i-l)  *9 
®  1  mY  <  J ,  4 )  — y  ( j ,  1 ) 

•2-del <2*nl < j,4) )-del <2*nl < j, 1) > 

••-•2/el 

print  t, "strain  at  tip",ee 
i-f  («•. gt.  dcri t)  be  < i)  -1 
40  continue 
return 
•nd 


subroutine  inputKhh.w) 
dimension  hh(10),w(20) 
read  (*,  *)  (hh  <i ) , i-1, 10> 
read  (*,  *)  (w(i ) ,  i«l,20) 
return 
end 


subroutine  gennod (x,y,nl,hh, w> 

dimension  x (200, 4) , y (200, 4) , n 1 (200, 4) , w (20) ,hh( 10) 

integer  col 

it=0 

nc-1 

k*r 

do  500  col-1, 19 
i f  (i t.  ge.  12)  k-4 
if (it. ge. 16) k*5 
if  (it. ge.2l) k»6 
if  (it. ge. 27) k=7 
if  (it.ge.34)k**8 
if  (it.ge.42) k-9 
do  20  i  - 1 ,  lc 

— - it-it+l-._  -  -  -• 

l=k+l 

8 (it, 1) (col > 
x (it, 2) *w (col+1 ) 

8 (it,3)=x (it, 2) 

x  (it, 4)*»x  (it,  1 ) 

y  ( i  t ,  1 ) =hh ( i  1 ) 

y(it,2)=y(it,l  ) 

y  (it,3)=<hh  (i ) 

y(it,4)»y(it,3) 

nl (it, l)-nc+i 

nl (it,2)-nl (it, 1)+1 

nl  (it,3)-nl  <it,2)~l 

nl(it,4)-nl(ittl)—l _ _ _ _  ... 

20  continue 

nc-nc+1 

500  continue 
return 
end 


subroutine  f ormn (n,x,y, nu, e,h, it > 
real  n,nu 

dimension  n(8,  8),x  (200,4)  ,y  (200,4) 
do  5  i«l,B 
n(5,i)-0 
n(l,il«0 

n(i,5)«0  : 

n(i,l)«0 
continue 
cl-1 
c2-l 

t—e/ (1— nu**2) 

w«t*  <l-nu) /2 

hK-K  (it ,2) -x  (i t ,  1 ) 

hy-y  <it,4)-y(it,  1) 

hx2-x (it, 2) **2-x (it , 1 ) **2 

hy2-y(it,4)ss2-y(it, 1) *S2 

hK3-K  (it, 2)  «3-x  (i  t,  1)  x*Z 
hy3-y(it,4)*«2-y(it, 1) **3 
n(2,2)-(t*hx*hy*h) /(cl**2) 
n  <2.3)-0 


n (2, 4) -( <t/2>  *hx*hy2*h)  /  <cl**2> 
n  <2.  6) *0 

n  <2,  7) * (nutt *hx*hyth) / (cl *c2) 
n (2,  8) -( ( (nu*t) /2) *hx2*hy*h) / (cl*c2) 
n (3,  3)  - (wf  hx  *hy*h ) / (c2**2> 
n  (3, 4) - ( (w/2) *hx2*hy*h) / (c2**2) 
n  <3, 6) - <w*hx*hy*h> / (cl*c2> 
n(3,7>-0 

n (3,8>-((w/2)*hx*hy2*h>/(cl*c2> 

n (4,4>-( ( (t/3) *hx Xhy3) / (c 1**2) + ( (w/3> *hy*hx3> / (c2**2) ) *h 
n  <4,  6) - ( (w/2) *hy*hx2*h) / (cl *c2> 

n(4,7>-( ( (nu*t) /2) *hx*hy2Xh) /(cl*c2)  ..  .  - 

n  <4, 8) * ( ( < (nu*t) /4)+(w/4> > *hx2*hy2*h> / (cl *c2> 
n(6,6)-(w*hx»hy*h)/(ci**2> 
n(&,7)-0 

n (&. 8)  —  ( (w/2) *hx*hy2*h> / (elf *2) 
n (7, 7) -(t*hx*hy*h) / (c2*X2> 
n (7, 8) - ( (t/2) *hx2*hy*h> / (c2*«2) 

n (8, B) - ( ( (t/3) Xhy*hx3> / (c2**2)+ ( (w/3) *hx*hy3> / (c 1**2) >*h 
do  10  i “2,  8 
do  10  1*2, i-1 
10  n(i,l)-n(l,i) 

return 
end 


— subroutine  mmult (ct,n, c, s, it) 

TT  '•  ~  dimension  ct  (B,  8) ,  c  (8,  8)  ,  s  (200, 8, 8) 

-  —  reel  n (8, 8) 

integer  it 
do  20  i-1, 8 
do  20  i-1,  8 
s ( i t , i , j ) -0 
do  10  k-1, 8 
do  10  1-1,8  : 

s(it,i,  j)-«(it,i,  j)-*-ct  (i,  k)*n(k,l)»c(l,j) 

10  continue  •  - 

20  continue 

.  return  .  : — „ — 

end 
c 
c 

subroutine  sti-fmat  (k,  s,  nl ,  n2,  it,  nen,  nei ,  nd-f ,  nei  > 

dimension  s  (200,  B,  8) ,  nl  (200, 4 ) ,  n2  (4) ,  n3  (4)  _ _ 

-  real  k  (330,330)  -  J: 

-  integer  it,nen,nel 
-  -  -  nn-nenXnd-f 

do  3  i *1,4 

3  n3  (i  )  — nl  (it,  i ) 

do  20  i-1, nei 

n2(i)«nd-f*(i-l) 
n3  (i  )  -nd-f  t  (n3  (i )  -1  > 

20  continue 

do  30  i-1, nei 
n j»n3(i > 
nl  -n2(i  ) 
do  40  i-1, nei 
nk-n3( j) 
nm-n2<i > 
do  30  1-1,  nd-f 
do  20  m-l,ndf 

k(nj+l,nk-*-m)-k(ni+l,nk'*m)*s(it,nl-*l,  nji+m) 
30  continue 

40  continue 
30  continue 
return 


c 

c 


end 


subroutine  boundcon < k , di sp , p , bw, be ) 
dimension  p (350) 
reel  k <350,350) 
integer  be (11) 

:«99  e  8 
do  10  i*l,350 
p<i>*0 
10  continue 
k  < 1 , 1 )  «z 
k<3,3)-z 
k  <5.5)*z 
k<7,7)»z 
k (9,9) *z 
k ( 17, 17)“z 
k  <  25 , 25 ) *z 
k < uu , 3u ) *z 
k  <43, 43) *z 
k<55,55)~z 
k  <69, 69) *z 
k  <85,85) ®z 
k  <103, 103) *z 
k<123, 123) *z 
*  -  k <143, 143) *Z 

k ( 163, 163) -Z 
k  <103, 1 83 ) *z  ' 

k  <203, 203) *z 
k  <223, 223) =z 
k  <243, 243) =z 
k  <263. 263)  =z 
l<  <283,283)  “Z 
k  <303, 303) *z 

k  <303, 305) =z  . . 

k <307, 307) «z 
. .  _  k  <309,309)  -z 

k  <31 1 , 3 1 1 )  »z  - - 

k  <3 13,313)  “Z 
k  <315,315)  *z 

k  <317,317)  mz  ‘ 

k  <319. 319)  *z 
k  <321 , 321 ) *z 
p  <2)*disp  -  ..  - 

cell  elim<k,2) 
p  <10)*diso 

call  elim<k,10)  -  .  ... 

p<18)*diso 

call  elim(k, 18) 

p  <26 >-disp 

call  elioCk, 26) 

p<34)»disp 

call  elim<k,34> 

p<44)«disp 

call  el im <k, 44) 

p  <56)*disp 

call  el im <k, 56) 

p  <70)«disp 

call  elim<k,70) 

p  <B6) *diso 

call  el im (k, 86) 

p  <104)*dlsp 

call  elim<k,104) 

p  <124>»disp 

call  elim(k,l24) 

0  <144)-diso 


cal  1  elim(k,  144  3 
?  < 164) «di sp 
call  elim<k,164) 
p  < 1B4> «disp 
call  cl im (k, 184) 
p  <204>«disp 
call  elim(k,204) 
p  <224) -disp 
call  cl i m  <k, 224) 
p  <244)«disp 
call  cl i m (k ,244) 
p <264) -disp 
call  dim (k, 264) 
p  < 284)— disp 
call  elim(k,2B4) 
p  <304) —disp 
call  elim<k,304> 
do  20  i-61, 161, 10 
j-int  <i/10)-3 
i-Mbc<j).lt.l)k<2*i,2»i)-= 

20  continue 
return 
end 

c„'  ‘  '  "  . 

c 

subroutine  elim<k,a) 
integer  a 
real  k (350, 350) 
do  10  i«l ,350 
10  k<a,i>«0 

k(a,a)-1.0 
return 
end 
c 
c 

subroutine  matsol v  <k, nen,  p ,  nd-f ,  del  > 
integer  nen,ndf 

real  k <330, 330)  - . 

dimension  del <330) , p (330) , indx <350) 
nn— nd-f*nen 

call  ludcmp <k, nn, 330, indx , q) 
call  lubksb (k,nn, 330, indx, p ) 
do  10  i-l,nn  . 

10  del  <i  )-p  <i  ) 

return  — 

end  ‘  • 

c 

c 

c 

subroutine  invert <a, n,np, indx, b,d> 
dimension  a(B,S),z(B,B) , indx  (B) 
do  12  i-l,n 
do  11  2»l,n 
c<i,j)-0 

11  continue 
=<i,i)-l. 

12  continue 

call  ludcmo<a,n,np,indx,d) 
do  13  j*i,n 

call  lubksb <a, n,np, indx , z < 1 , j > ) 

13  continue 

do  14  i-l,n 
do  14  j-l,n 

14  a<i,  j)-s<i,  j) 
return 

end 


* 


11 


13 

14 


13 

16 

17 


ia 

19 


subroutine  ludcmp (a, n , np , indx , d ) 
dimensi on  * (np , np ) , i ndx (n) ,  vv (330) 
tiny*l.O  e-10 
d-l. 

do  12  i*l , n 
aamax “O . 
do  11  j«l,n 

if(abs(a(i,j>>.  gt. aamax )  aamax “abs  (a(l,j)  ) 
continue  '  ' 
if (aamax . eq. 0. ) aamax*tiny 
vv (i ) «l . /aamax 
continue 
do  19  j*l.n 

if ( j.gt. llthen 
do  14  i«i,j-t 
sum=a (i , j ) 
if (i .gt. llthen 
do  13  k«l,i-l 

sum*sum-a ( i  ,  k )  *a  (k  ,  j  > 
conti nue 
a (i , j 1 “sum 
endif 

■  —  continue- - - 

endif 


a*max**0. - *-  - - - 

do  16  i*j,n 
sum*»a(i  ,  j) 
if  (j.gt. 1)  then 
do  13  k*i,j— 1 

sum“sum-a (i,k)*a(k,j) 
continue 
a ( i , j ) =sum 
endif 

dum=w(i  )  Sabs  (sum) 
if (dum.ge. aamax) then  ' 
imax“i 
aamax<*dum 
endif 
continue 

if  < j . ne. imax ) then 
do  17  k«l,n 

dum=a(imax, k) 
a  (imax,  k)“a  ( j  ,  k> 
a< j, k)»du« 
continue 
d“-d 

w  ( i  max ) *  w  ( j ) 
endif 

indx tj)“i max 
if  (j'.ne.n)  then 

if (a( j,  j)  .eq.O. >a(j, j)“tiny 
dum»l.  /a  ( j-,  j) 
do  18  i“j*l,n 

a  (i  ,  j )  «a  (i  ,  j)  *dum 
conti nue 
endif 
continue 

if  <a(n,n).eq.0)a (n, n) “tiny 

return 

end 


c 

c 


subroutine  lubksb (a,n,np, indx.b) 
dimeneion  alno.no) .indx (n> .b <n> 


11 


ii*0 

do  12  i«i,n 
1 l»indx <i > 
sum»b (1 1) 
b(ll)-b<i) 
i  -f  C i  i .  ne.  0) then 
do  11  j*ii,i-l 

sum=sum-a (i ,  j  )  *b  ( j  ) 
conti nue 

els®  i f (sum. ne. 0. )  then 
i  i  =i 
endi-f 
b (i )«sufli 

12  continue 

da  14  i«n,l,-l 
sum*b ( i ) 
if  (i  .  lt.n)  then 
do  13  j^i+ljn 

sum»sum-a (i,j)fb(j) 

13  continue  . 

endi-f 

b  Ci  )  =*sum/a  (i  ,  i  ) 

14  continue 

return - —  - - 

®nd  ...  . .... 


subroutine  -f  ormcinv  (c . s« , y , et , it , ck) 
dimension  c (B, B) ,x (200,4) , y (200,4) ,ct (B,B) 
dimension  temp  (8,8),ck(5,8,8),  indx  (8) ,  b  (8) 
do  2  i«l,B 
do  2  j-1,8 

2  temp ( i , j ) =0 
j=0 

do  5  i*l,8,2 
j*j+l 

temp (i, 1)«1 

temp  (i  ,2)«x  (it,  j  )  —  -  ;.  _ 

temp (i,3)»y(it,  j)  •  -  - 

temp(i,4)-x(it,  j)*y(it,j) 

3  continue 

j-0 

do  10  i-2,8,2 

j-j+1 

temp (i,3)*l 

tempi!  ,6)»x  (it,  j)  -  - 

temp  (i,7)»y(it,  j) 

temp  (i  ,B)»x  (it,  j)*y  (it,  j) 

10  continue 

call  invert (temp, 8,8, indx,b,d)  .  . 
do  20  i»l,B 
do  20  j-1,8 

20  c  (i ,  j)» temp  (i  ,  j) 

do  SO  i-1,8 
do  30  1—1,8 

ct ( i , 1 ) *c ( 1 , i ) 

i-f  (it.  ge.37)  ck (i t-36,  i  ,  1  )-c <i  ,1  > 

SO  continue 
return 
end 


subroutine  graphics (x, y, ni ,nel , del , tx ) 
dimension  x (200,4) ,y(200,4) ,nl (200, 4), del (330) 
Inteqer  xo.yol.vo2.xo2.vo3.vc4 


„  do  20  i*l,nel 
do  10  j«l,4 
KO“i  j: 
yol*iyl 
yo2*iy2 
::t«0 
yt“0 

t*x  <i ,  j )  +del  <2*111  <i  ,  j > ~1 ) 
yt*y (i , j ) +del (2*nl (i , j ) ) 

<>it*30+300) 
iyl-int (400-yt*50> 
iy2*int  <400-ryt*50) 
call  plot  (i:s,iyl> 
call  plct(ix,iy2> 
if (j.eq. Dgo  to  3 
call  draw(xo,yol, ix, iyl) 
call  draw <xo, yo2, ix , 1 y2) 
i-f  ( j  .  eq .  4)  go  to  7 
go  to  10 
3  ::o2=iv. 

yo3=iy 1 
yo4-iy2 
go  to  lO 

_ 7  ....  call  draw (ix, iyl , xo2, yo3)  -  - - 

call  draw (i x , i y2, xa2krn ( ) 

return  — . 

end  -  •  ~ .  r—  ...  •:  -• 

c 

c 

subroutine  forme <del . e, nu, n  1 ,  y , nec, it > 
dimensi on  del (350) , n 1 (200, 4) , y (200, 4  > 
real  lm, 1 , nu, nec (200) , Ub 
t»300 

kb»1.38  e-23 

g=nec (it ) *kb*t  -  ■  r' 

i=y(it,4>-y (it, 1) 

lm-(l+del <2*nl (it, 4) )-del <2*nl (it; 1) > > /I  - 

•f  l="func<lm)  — - - *  -- 

f l*f l*g*(lm-(l/lm)**2> 

lm»(  (1-rdel  (2*nl  (it,  4) ) -del  (2*ni (it,  1)  >  )/l  1+0.2 
f2-func(lm> 

f2*f2*g*(lm~(l/lm>**2) 

•«(f2~fl)*5 

return  r  .  - 

end 
c 
c 

■function  -funcllm) 

real  In, lml 

lml»l/lm 

n«100 

a«1.4 

b-4.B4B98  a  -2 
e«5. 38367  •  -3 
C-2.70111  e  -4 

f  l*l«-a*<l/n)*<3*lm*x2+4*li»>l> 
•F2«b*<a/n)**2)*(3*lm**4*B*lm+B*lmI**2>  ■ 

f  3»c *  ( <  1  /n )  * *3)  *  <33*  1  m**6+60*  1  m* *3+72+64*  1  m  1  f  *3) 

■?4«d*  ( <l/n>**4)*  (630*1  m**B+ 1 120*1  m**!^  1440*1  m**2*  1336*1  ml 
1  *1280*iml**4) 

•func*f  l-*-f  2-rf  3-*-f4 
return 


end 


subroutine  calcnec  (nec, del , dt , nl , y, it) 
real  nec (200) ,del (350)  ,  dt,y (200,4)  ,1 
integer  nl <200, 4), it 
1-y  (it,4)-y  (it,  1) 

dl«del <2*nl (it, 4) >-del (2*nl  (it,  1) > 

stn®dl /I 

al-f-1 

nec  < it ) =nec (it)  *exp  (-al-f  *stn*dt  > 

return 

end 


^include  <stdio.h> 

#include  <sys/gra-f .  h> 
main  C) 

{ 

struct  Form  screenForm;  - 

i-f  <g_init  (S<screenForm)  ) 

<  print-f  ( "  ERROR*  G_init  -failed\n">5 

_  y 

g_cl  screen  (>  f 

-  g_-flush(>;  — — 

>  *  .  . 

plot_(x,y) 
int  *:<,*>•} 

g_writepixel (»x, *y, Black) ; 

"X 

* 

newskrn  (  ) 

{ 

g_-flush()j 

% 

a  •  —  • 

draw_ (xO, yO, x 1 , y 1 ) 

int  *x0,*y0,*xl,*yl5-  - -  - 

g_drawl ine (4x0,  *y0,  *>:  1 ,  *yl )  ; 


10,0. 1, 1.0, 0.5,0. 1, 132,4,2,0. 1,1. 0,0. 8, 0.5, 4 
6.0,4.0,u.0,2. 0,1.B,1.6,1.2,0.  8,0.4,0.0 
-6. 0,-4. 0,-3. 0,-2. 0,-1. 5, -1.0, -0.6, -0.4, -0.2 
0.0, 0.2, 0.4, 0.6, 0.8, 1.0, 1.5,2. 0, 3. 0, 4 . 0, 6. 0 


